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EDITORIAL 


Beginning with this issue the Proceedings will appear in two Parts, 
“A” and “B”. This change has been dictated by the increasing 
volume of important communications which has proved too great even 
for the larger amount of space available since monthly publication began 
in January 1948. 


It is intended provisionally to publish in Part A papers dealing 
with subjects such as_ the quantum theory, statistical. mechanics, 
nuclear physics and cosmic rays, atomic physics, molecules, spectra, 
theories of solids, liquids and gases, surface phenomena, growth and 
properties of crystals, crystal structure, luminescence, electrodynamics, 
heat and thermodynamics, and standards. Part B will include subjects 
such as acoustics (including ultrasonics), optical design, electron 
optics, colour, elasticity and other mechanical properties of solids and 
liquids, crystal structure analysis, magnetic materials, refrigeration and 
liquefaction, electric discharges, radio, geo- and ionosphere physics, 
astrophysics, and solar physics. It is obvious that a hard and fast 
dividing line cannot at present be drawn, but we feel sure that in a short 
time each Part will develop a character of its own, so that the reader 
will have little difficulty in deciding with which he wishes to be supplied. 


With this new arrangement we hope to be able, not only to give 
improved service to Fellows and other subscribers, but also to effect 
that speedy publication of important scientific papers which is so vital 
to the development of physics. 

G. INGLE FINCH, 
President. 
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On the Theory of Aplanatic Aspheric Systems 
By G. D. WASSERMANN *#® anp E. WOLF t 


H. H. Wills Physical Laboratory, University of Bristol 


* Now at King’s College, University of Durham 
+ Now at the Observatories, University of Cambridge 


MS. received 21st February 1948; read before Optical Group 17th December 1948 


ABSTRACT. Methods are described for the design of two aspheric surfaces for any given 
centred system, so as to achieve exact aplanatism. The practical application of the methods 
is illustrated by the design of a reflecting microscope. 


St) IND ROD ELTON 


HE design of aspherical systems raises a number of questions. Foremost 

among these is the problem of how to determine for any given system the 

aspheric profiles so as to obtain good definition over the whole of a finite 
field. To the approximation of third order (Seidel) optics this can be accomplished 
by the elegant method of plate diagram analysis (Burch 1942, Linfoot 1944). As 
a rule, however, higher aberrations must be taken into account. A reasonable way 
of doing this is to supplement the preliminary Seidel design by further analysis. 
G<neral methods using differential corrections were described by Volosov (1947), 
but their form makes them difficult to handle. Frequently, however, it is 
sufficient to confine the further analysis to the attainment of either axial stig- 
matism or aplanatism.{ The former has been discussed by Herzberger and 
Hoadley (1946) and by Wolf (1948). No corresponding analysis for the case of 
aplanatism appears to have been made. However, a simple method adequate for 
some systems was given by Linfoot (1943). Particular systems were considered 
by various authors; for instance, Schwarzschild (1905) and Chrétien (1922) 
deriyed exact equations for two-mirror aplanatic telescopes. Their methods, 
although valuable for the practical design of such aplanats, are not very 
suitable for application to more complex systems.§ 

The aim of the present paper is to provide methods for the design of two 
aspheric surfaces for any given centred system, so as to achieve exact aplanatism. 
We derive two simultaneous first-order differential equations for their profiles. 
‘These equations contain certain functions which characterize two-ray congruences 
in the spaces immediately before and after the two surfaces. These are easily 
obtained from a ray trace. Our differential equations permit of a rigorous 
solution of the problem and can be solved to any desired accuracy by standard 
numerical methods. 


§2. EQUATIONS FOR THE DESIGN OF TWO SURFACES 
TO ENSURE APLANATISM 
2.1. In this section we derive formulae for the design of two correcting 
surfaces so as to obtain simultaneously: (a) axial stigmatism, (b) exact satisfaction 
of the sine condition. Our results apply to any centred system in which the two 


¢ We use the term aplanatism to imply axial stigmatism together with the satisfaction of the 
exact sine condition. 
s § This is illustrated, for exarnple, by the analysis of Bureau and Swings (1934) who applied 
‘Sc}:warzschild’s methods to the design of two-mirror aplanats with both foci finite. 
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! aspheric surfaces are optical neighbours. The two surfaces in question may, 
_ however, be separated from the object or image points by any number of spherical or 


aspherical surfaces. We shall only be concerned with the final corrections of the 
system, and we assume all design data, save the profiles of the two correctors, to be 
known. : 

We introduce two sets of Cartesian coordinate systems (Figure 1) with origins at 
the poles O and O’ (assumed to be known) of the two correcting surfaces & and 2’, 
with their x axes along the axis of the system; and take OO’ =d. Let us denote by 
o and o’ the spaces immediately preceding » and following D’, and by o* the space 


- between them. F inally, let n, n* and n’ be the corresponding refractive indices. 


Figure 1. 


‘The rays which proceed from the axial object point S are refracted (or reflected) 
at the successive surfaces and form normal rectilinear congruence I in the space o. 
Let h denote the height at which a typical ray of I meets OY (if produced far 
enough) and let w denote the angle which the rays make with the w axis (measured 
anticlockwise from the positive x direction to the direction along which the light 
advances). I’ can then be completely specified by two relations: 


w=a(t); PPSTEEN, ae SAP, EL ERE Nee se (2.1) 


where ¢ is a free parameter.’ If the axial object point S is at a finite distance we 
choose 


BE SUINOT apie hPa na oe (2.2) 


where © is the angle which the corresponding ray in the object space makes 
with the x axis (Figure 1). If, on the other hand, the object is at infinity, we 


choose t=H where A is the distance of the corresponding ray in the object 


space from the x axis (Figure 2(a)). 


Figure 2 (a). Figure 2 (0b). 


The ray specified by the parameter ¢ will be referred to as the ray ¢ and its 
intersection with X will be denoted by T. As the relations (2.1) can rarely, in 
practice, be obtained in analytical form, we form a table of w and h based on a ray 
trace for different values of ¢. Some designers would probably favour poly- 
nomial approximations to w(t)and (2) at this stage ; but the merits of this expedient 
are doubtful and will be discussed later. 

1-2 
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Similarly the rays which are traced backwards from the axial image point S° 
are refracted (or reflected) in succession by all surfaces which follow ©’ and form a — 
normal rectilinear congruence in the space o’. Let I” be the congruence in o’ 
which is obtained by reversing the sense of every ray of the “back-traced” 
congruence, i.e. we regard the rays of I” as though they proceed from &’ towards 
the axialimage point (Figure 1). I’ can then be specified by two relations: 


Weal. hah ne oe (2.3) 


where w’ and h’ are now referred to axes at O’ and the parameter ?¢’ has a signifi- 
cance analogous to that of t (Figures 1 and 2(b)). The quantities w’ and h’ are — 
obtained for selected values of t’ by tracing rays backwards from the axial image 
point S’ into the space o’. . 

The conditions (a) and (d) for aplanatism require that the ray ¢ of I be trans- 
formed by two successive refractions at & and =’ into the ray ¢’ of I’ where t and ¢” 
are related by the sine condition, which takes the form 


t/t = CODStANE 5 pe (2.4) 


Equation (2.4) defines uniquely the ray t’ corresponding to a given ray ¢, so that our 
problem can be regarded as involving only the free parameter ¢. To avoid 
cumbersome suffixes we shall use undashed and dashed symbols to refer to quan- 
tities associated with a pair of corresponding rays. 

2.2. Let («, y) and (x’, y’) be the coordinates, referred to axes at O and O’ res- 
pectively, of the points T and 'T’ (Figure 1). We have, by Snell’s law, for the 
refraction at & 


n(cos w dx/dt +sin w dy/dt) =n*(cos w*dx/dt+sinw*dy/dt), ...... (2.5) 


where w* is the angle which TT’ makes with the axis. 
We have from Figure 1 


sinwo* = KojR; - cos wie Re | eee (2.6): 
where : 
R, =x —x-+d,» R,=y'—y, and R= RAR 2. oss (2.7) 
Also, from the figure, 
veh -extaniog ae Si eee a eee (2.8) 


and y =h'+x' tana’. 


Substituting for cos w* and sin w* from (2.6) and for dy/dt from (2.8) into (2.5) it 
follows that 


. 
. 
| 
. 
dx (n* R,,—nR cos Lys dh | 
pore : tne | eee stan |. 


(n*R,,—nR sin oe 


. De it a cag ec SS a ieee ee (2.10) 
Similarly we obtain 
dx' ___[(n*R,—n'Reosw’) ae ; 
dis (n “Kyi hsina) pra qian |. 
Antes 249 


Equations (2.10) and (2.11) permit together with (2.7), (2.8) and (2.9) a complete. 
computation of Yand&’. For by means of (2.7), (2.8) and (2.9) we could eliminate. 
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y and y’ from (2.10) and (2.11) and thus obtain two simultaneous first-order 
differential equations of the type 


Gyidi=fir rot). and, Ox |Gi= p(n, eet)e 2 eee, (2:12) 


‘These may be integrated numerically by standard methods. But as we require y 
and y’ as well as x and «’ for a range of values of the parameter ¢, it is preferable 
to avoid elimination and solve for the unknown quantities step by step. - This will 
be discussed in more detail in the next section. Equation (2.12), subject to the 
_ boundary conditions x =x’ =0 for t=0 show that if there exists a physical solution 
of our problem it will be unique. 
2.3. Ifthe correcting surface ¥ is a mirror, the fellawing rule should be applied : 

Put n= —n* =1 and change R,, and R, in (2.10) into —R, and —R,. Similarly 
if X’ is a mirror, put n’ = —n* si and Shane the sign of R, and R,, in (2.11). 


§3. METHODS OF COMPUTATION AND SOME APPLICATIONS 


3.1. We shall first summarize the procedure for computing the two correcting 
surfaces. 

(a) Rays are traced from the axial object point into the space o for a set of 
selected values of the parameter ¢. For each value of ¢ the quantities w and h are 
determined. Following this we calculate the values of dtanw/dt and dh/dt by 
numerical differentiation or otherwise. In simple cases (e.g. when & is the first 
surface of the system) the differentiation may be performed analytically. 

(6) For each value of ¢ the corresponding value of t’ is calculated from (2.4). 
The quantities w’ and h’ are now determined by tracing rays backwards from the 
axial image point into the space o’. The differential coefficients dtan w’/dt and 
dh'/dt are subsequently determined as above. 

(c) Next, equations (2.10) and (2.11) are integrated numerically step by step, 
for every value of x and x’ the corresponding values of y and y’ are calculated at 
every step from (2.8) and (2.9). The integration is performed in two stages. 
‘The first few values can be obtained by the application of the methods of Runge 
and Kutta (Runge and Konig 1924, p. 286). The integration can then be con- 
tinued by Adam’s method (Whittaker and Robinson 1946, p. 363). In this way 
we obtain the two profiles in tabulated form. 

Some opticians would probably employ throughout the design polynomial 
approximations, based on the tracings of a small number of rays. Although 
polynomial approximations form an important tool in preliminary design they are, 
in our opinion, not necessarily of advantage in the final computation of an aspheric 
system. Our doubts of the expediency of polynomial approximations arise first of 
all in connection with the high accuracy which is frequently required in the design 
of an aspheric surface. (In astrographic cameras for example, it is generally 
desirable to determine the profiles down to a fraction of a fringe.) On account of 
slow convergence (cf. Martin 1944, p. 108) this may necessitate the evaluation of a 
large number of coefficients. Asa rule, however, it is difficult to obtain more than 
the first three or four coefficients in terms of the design parameters, so that the 
expansions have to be supplemented in particular cases. Moreover, if we base the 
expansion on the tracing of a small number of rays we encounter further difficulties 
as it is generally not possible to determine in advance the number of terms re- 
quired for a satisfactory expansion (Herzberger and Hoadley 1946, p. 339).. 
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In view of these considerations it is our opinion that polynomial approximations _ 
are not very suitable in the solution of this type of problem and should be replaced _ 
by tabulated functions based on tracings of dense fans of rays. Linfoot (1943, 
p. 494) also advocates this procedure. ; 

3.2. ‘To test the analysis we epplied our methods to the design of an apJanatic 
reflecting microscope and compared the 
results with those obtained from the 
analytical solutions of Schwarzschild (1905) 
and Chrétien (1922) for two-mirror aplanats. 
Microscopes of this type were constructed 
and discussed by Burch (1947). _ 

The objective (Figure 3) consists of a 
concave primary and a convex secondary 
mirror. ‘These are separated by a distance 
d and image parallel light at a distance 
m from the pole of the primary (mandd are negative in Figure 3). In this case the 
ray traces are very simple. In the previous notation we have 


Figure 3. 


w=0, h=—mtanO; wy 0), = hee es (3.1) 
Corresponding parameters satisfy the sine condition which now takes the form 
Biz fh ang eel) on EEN eae eee (3.2) 
where 
t=sinQ@; Pas He Pee oe (3.3) 


and f is the focal length of the system. On applying the rule for mirrors of 
§2.3 and making use of the above relations, equations (2.10) and (2.11) take the 
form 


d cose = 
= =— i ooo = + tan | (x= m)\secO <7. (3.4). 
ax [die fRA(R= RR). Pia oo, 8s eens (3.5) 


where as before 
R,=x —x+d, R,=y'—y, and R®*=R,*+R,?. 


PRS (3.6) 
Also (2.8) and (2.9) become 
y=(x—m) tan@- 97> © 3 Paes (3.7) 
and 
y’ =fsin Cree i NG hs (3.8) 


We computed the system with design data m= —6cm., d= —3cm., f=0-75cm. 
and numerical aperture (N.A.) =0-58. 


Twenty equally spaced values of the parameter ¢ were selected, viz. 
ma O0297 eal oe). 


First the values of x, and x,’ with the corresponding values of y, and y,’ were 
computed for r=1 and r=2 by Runge’s method. The integration was then 
continued by Adam’s method. In calculating the first few differences occurring 
in Adam’s formulae use was made of the fact that Mi axsand H.p=%./. 0 Results 
for the first six values of x are shown in table 2, where A denotes the difference 
between these solutions and those obtained from Schwarzschild’s formulae. 
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Owing to the inaccuracies of Runge’s method an almost constant error of 
— 8 x10-%cm. occurred in all values of x, for 7>2. The maximum error in x,’ was 
2 x10-*cm. 


Table 1. Reflecting Microscope, m= —6cm., d= —3cm., f=0-75cm., 


N.A.=058. 

Concave primary (in cm.) Convex secondary (in cm.) 
rT t=sin ©), =X, OSAoe. Vp 108Ay,. en HOS ame eee 
ofl 0-029 0-003469 0 0-173973 0 0-000276 0 0:02175 
2 0-058 0-013880 0 0-347781 0 0-001104 1 0:04350: 
3 0-087 0-031240 0 0:521258 1 0-002486 0 .0:06525 
4 0-116 0-055562 1 0694241 1 0-004424 0 0:08700: 
5 0-145 0:086866 1 0:866562 1 0-006920 0 0-10875 
6 0-174 0-125175 0 1-038054 0 0:009978 1 0-13050: 
ih 0-203 0-170523 0 1:208547 0 0-013603 1 0-15225. 
8 0-232 0-222949 1 1-377870 0 0-:017799 —1 0-17400 
9 0-261 0-282499 1 1-545849 0 0-022575 —1 0-19575 

10 0-290 0:349231 0 1-712307 0 0-027937 —1 0:21750 
11 0-319 0-423213 0 1:877063 0 0-033894 0 0-23925 
12 0-348 0504524 0 2-039933 0 0-040456 0 026100: 
13 0-377 0:593259 0 2200726 0 0-047634 0 0-28275 
14 0-406 0-689530 0 2°359245 0 0-:055441 0 0-30450 
15 0-435 0:793468 —1 2°515288 —1 0-063892 1 0-32625 
16 0-464 0-905230 0 2668638  —1 0-073002 2. .0-34800 
ily 0-493 1-024999 —1 2:819072 0 0-082790 1 0:36975 
18 0-522 1:152996 aw 2:966351 0 0:093277 —1 0-39150 
iLO 0-551 1-289480 =) 3°110222 —1 0-104487 0 0-41325 
20 0-580 1-434771 —1 3250403 1 0-116445 —{ 0-43500 


A denotes difference between results obtained from Schwarzchild’s formulae and those 
calculated by the present methods. 


able 2 
7 1 2 3 4 5 6 
Se 0-003471 0-:013884 0-031248 0-055570 0-086874 0-125183 
108Ax, 2 4 + 8 y 9 8 


To obtain higher accuracy we repeated the integration for the first two points. 
This was carried out by substituting the values of dx,/dt and dx,'/dt (r=0, 1, 2) 
previously obtained into the formulae of Newton-Cotes (Whittaker and Robinson 
1946, p. 152). Very close agreement with Schwarzschild’s results was obtained. 
The integration was then continued as before with results givenin Table 1. It can 
be seen that A never exceeds 2 x10-®cm., i.e. about 1/20 of a reflection fringe. 
Such high accuracy in the design data is not needed by the practical optician, 
but it is a necessary preliminary to a numerical analysis of the off-axis errors, 
on which any theoretical assessment of the value of the design must be based. 

The computation was carried out with the aid of Gifford’s 8-figure tables 
of trigonometric functions (interpolating up to 1/100’’) and Milne-Thomson’s 
8-figure tables of square roots. 

Since, in the derivation of our equations, no use was made of any approxi- 
mations, any inaccuracies in the final solution are due to small computational 
errors. 


8 W. Ehrenberg and R. E. Siday 


ACKNOWLEDGMENT 
We are indebted to Drs. C. R. Burch, H. Frohlich and E. H. Linfoot for 


helpful criticisms and encouragement. 


REFERENCES 


Burcu, C. R., 1942, Mon. Not. R. Astr. Soc., 102, 159; 1947, Proc. Phys. Soc., 59, 41. 

Bureau, F., and Swincs, P., 1934, Rev. d’Optique, 13, 127. 

CuriTIEN, H., 1922, Rev. d’Optique, 1, 49. 

HeERzBERGER, M., and Hoap ey, H. O., 1946, 7. Opt. Soc. Amer., 36, 334. 

Linroot, E. H., 1943, Proc. Phys. Soc., 55, 481; 1944, Mon. Not. R. Astr. Soc., 104, 48. 

Martin, L. C., 1944, Proc. Phys. Soc., 56, 104. 

Runcg, C., and Konic, H., 1924, Numerisches Rechnen (Berlin: Springer). 

SCHWARZSCHILD, K., 1905, Géttingen Observatory. 

Votosov, D. S., 1947, ¥. Opt. Soc. Amer., 37, 342. 

WuitTakER, E. T., and Rosinson, G., 1946, The Calculus of Observations (London : 
Blackie & Son). 

Wo tr, E., 1948, Proc. Phys. Soc., 61, 494. 


~ 


The Refractive Index in Electron Optics and the 
Principles of Dynamics 


By W. EHRENBERG* anp R. E. SIDAY tT 
* Birkbeck College, University of London + I.C.1. Fellow, Edinburgh University 


MS. received 23rd February 1948, and in amended form 28th Fuly 1948; 
read 3rd December 1948 


ABSTRACT. In view of mis-statements made in the literature, the origin of the refractive 
index in electron optics is discussed in some detail, and the uniqueness of an expression 
previously given is demonstrated. On this basis, some general properties of electron 
optics are investigated. 

A relation between ray direction and wave normal is obtained. Whereas the refractive 
index is unique in terms of the magnetic vector potential A, this itself is arbitrary to some 
extent. It is shown that A must, for purposes of electron optics, be chosen so as to satisfy 
Stokes’ theorem and that, if it does, no observable effects result from the arbitrariness of A. 
An expression for the optical path difference is given in terms of the magnetic flux enclosed. 
The results are applied to a number of questions, viz. the differential equations for tra- 
jectories, the focusing properties of an axially symmetric field and the interference pattern 
produced by two converging bundles of rays which enclose a magnetic flux. 


wi LN TRODUCTION 


HEREAS in light optics the refractive index of a medium is in the first 
instance an experimental datum, and its accurate value is the basis of 
any detailed discussion of the performance of optical instruments, all 

geometrical electron optics is entirely contained in Lorentz’s equation for the 
forces acting on a moving charge. As a result, all equations for trajectories can 
be derived directly from that equation by specifying the electric and magnetic fields. 
Thus the réle of the refractive index in electron optics is far less obvious than that 
of its counterpart for light, and in fact different authors have proposed essentially 
different values of the refractive index for the same field without arousing much 
perturbation (Glaser 1933, 1937, Opatowski 1943). Only through the persistent 


ee 


Refractive Index in Electron Optics 9 


use of the refractive index, however, is the optical character of electron optics 
fully brought out, the flow of reasoning maintained and in fact greatly simplified, 
and the way opened for a wave optical treatment of electronic problems, such as 
the resolving power of the electron microscope. 

It is Glaser’s great merit to have first pointed out the existence of this refractive 
index and to have given its value correctly. The reasoning that has led him to it is, 
however, as will be shown, inexact, so that a logical extension of his line of thought 
has led to wrong values (Opatowski 1943). It is the object of this paper to review 
the source of the refractive index in electron optics and to discuss some of its 
general properties which in view of its peculiar indefiniteness require examination. 


§2. THE DERIVATION OF THE REFRACTIVE INDEX 


The term ‘electron optics” is often justified by reference to the corre- 
spondenc > between instruments designed for electrons and for light. A systematic 
definition of optics in a generalized sense appears to arise from ‘‘Fermat’s 
principle”. Although in its historical form it is a principle of fastest arrival, it is 
now recognized as 

Pip 

Colts eae meet el ly a See. (1) 

Pp, J 
Here ds stands for the length of an element of ray, ». denotes the refractive index, in 
general a function of position, direction and frequency, | wds is the line integral 
over any line between the limits indicated as subscripts and superscripts of 5.* 

Fermat’s principle may be enunciated as follows: In any optical medium a 
scalar quantity py, finite everywhere in space, can be defined so that the line integral 
in the three-dimensional space taken between any fixed points is an extremum for 
the ray which passes through these points. Or, more precisely: If the value of 


is pds along a particular line joining any two points Py and P, of an optical system 


is stationary with respect to weak variations that leave the end points unchanged, 
then the line is a ray of the optical system. Asan obvious corollary it follows that : 
If no line can be drawn joining P, and P,, such that fds taken along this line is 
stationary for weak variations, no ray can pass according to geometrical optics from 
Boito-P,. 

In the optics of light the value of » as a measured quantity 1s finite and single 
valued and is continuous except at a finite number of surfaces separating different 
media. As will become apparent later, however, more general mathematical 
expressions can be used in place of the simple values usually given, and in electron 
optics they arise in fact automatically. Particular expressions are admissible only 
if they do not violate the conditions implied for the validity of Fermat’s result. 
Firstly, the difference of Fermat’s integrals along any two paths must be defined 
anywhere in space where the optical properties are to be investigated. This 
requires that the refractive index should be fixed everywhere in space once it is 


* Sometimes the variational operator is written 6 for a variation in which the limits of the 
variable of integration are to be kept constant, and A if these limits should be varied. Sometimes the 
boundary conditions are not given in the equation but must be found from the context. We 
find this confusing and propose to write the complete boundary conditions at the lower and upper 
limit of the integral below and above the 5 sign, so that, for instance, Xo, Vo, %o; to below 5 denote these 
values for the line at the lower limit etc. x9, Vo, 20) to-+Ato denote that at the lower limit 
X=2X, V=V etc. but that the value for tp should be varied; or in other words, that the lines under 
comparison must go through xp, 9, %o, but may do so at any value t near ty. We shall also use the 
symbol Py to denote the set of coordinates x9, Vo, Zo and P, for x, V1, 21. ' 
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fixed in the neighbourhood of one point, so that » must be single valued. Secondly, 
it should have no singularities through which Fermat’s integral would become 
infinite. And thirdly, any discontinuities should be of such a nature that they 
appear as limiting cases of acontinuous p. This is true for the boundary between 
two optical media but would not be true if, for example, 1 were proportional to a 
cyclic coordinate which is limited to the values 0 to 27 in order to make it single- 
valued, for the numerical values of the integrals along either side of an irremovable 
discontinuity will not be equal even when the lines are infinitely close together. 
That would make the weak variation required by Fermat’s principle impossible.. 

Fermat’s principle does not contain time; and, in fact, not only in geometrical 
optics (Gabor 1945) but also in wave optics the timing of the phenomena observed 
is of no interest. Problems into which timing enters, such as the velocity of light, 
stand quite apart. ‘Thus Fermat’s principle suggests a generalized optics when- 
ever problems arise in any field in which timing is of no interest and a function 
of position and direction can be defined so that rays are determined by an expression 
of the form ot Fermat’s principle. 

It is now important to see that such a generalization cannot be based on 
Hamilton’s principle 

Betiwe 
Neo! es ee ee eS (2) 
Posty 2 

in spite of the fact that its integrand is readily transformed into the form of Fermat’s 
integrand. Hamilton’s principle may be enunciated: The world line of a particle 
(i.e. the line in the x, y, z, space) which goes through the points Po, ¢) and P), t, is 


“ty 
the line going through these points which minimizes the line integral | L dt. 
ae: 


its general validity arises from its equivalence with the equations of motion 
of Lagrange : ) 


Obs ds for 
ope alsa) aoa | Sra al bias eae (3) 
For electrons in particular, 
__ mv? e 
L= aes e® + 7 (7. 0 aaa sae Aes (4) 


(Schwarzschild 1903). It is only necessary to write this L in, say, Cartesian 
coordinates x, y, 2, x, y, 2, to perform the differentiations (3) and to insert 
—grad ®=E, curlA=H in order to arrive at Lorentz’ equation for the acceleration 
of an electron 


mv = eE + <[v xH] re (5) 


in Gauss’ units, with e= —4:°8 x107-!%£.s.v. 

That Fermat’s principle cannot be associated with Hamilton’s principle was, 
of course, clearly realized by L. de Broglie, but later authors who have specialized 
in electron optics have tried to associate the two principles by equating u = L/v. 
Only for a particular—valid but unjustified—choice of L does this lead to a correct 
answer. Glaser, in fact, by this artifice obtained a correct value, whereas Opatow- 
ski, by generalizing Glaser’s choice of L, was led to a set of values for 4 some of 
which are wrong. Glaser’s correct result appears to be due to such a combination 
of an arbitrary breaking up of his integrand into two terms with a’neglect of 
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Hamilton’s boundary conditions that Lagrange’s form of the principle of least- 
action (see below) is in fact substituted for Hamilton’s principle. In some text-. 
books on electron optics (e.g. Maloff and Epstein 1938) the correct variational. 
principle is used while in others it is not. We have found no explicit discussion 
of this point in any, and the choice of a principle seems to have been based rather on 
“Tuck than judgment’’. 

It is quite correct to transform the integrand Ldt since a given line x=</(t),. 
y=y(t), 3=2(t) which satisfies x) =x(to), Vo =(to), --- 2, =2(t,) and which gives 


bh 
to the integral | Ldt the value B, will give the same value to the integrals. 
Ji 


cP i ma 
| (L/v) ds or | (L/x)dx etc. so that Hamilton’s principle can alternatively be. 
Po v Xo 
written 
Pit) fe 
5 |(Lle)ds=0, oe (6) 


Po; bo 


th. Yrs 21, @1 
3 | (L/a)de=0 ete, ee (6a) 
ton don cortore 
where in (6a) L/x is now a function of x, y, 2, dy/dx, dz/dx, dt/dx. Further, it is. 
easily seen that equations (3) can be combined to 


(d/dt)(L—Zg(OL/0g,))-OL/fot=0- ae (7) 
so that 
ie MOI OG Const irs ie ics (7 a) 


(—h=the energy of the electron) if Z does not depend explicitly on ¢, i.e. if the- 
magnetic and the electric fields are constant in time—the case we are here: 
exclusively interested in. As 


R= ani ==ma/(l+y?-s*), vaegx = 2x, <y =dy/deete.,. cn. (8) 


we can eliminate the variables «, y, &, v from L/v by means of (7a) and (8) and_ 
replace them by y’, 2’, h so that the integral in (6) takes the form 


CDSE Oya nee peaelln ldnciae (9) 


This does not contain ¢ or derivatives of the form d/dt and is, in fact, a line integral’ 
in three dimensional space; that is, it depends only on the projection of the world 
line on three dimensional space, with / as a parameter, as the analogy to Fermat’s. 
principle would require. 

But the value of h in this integrand can be determined numerically only after’ 
Lagrange’s equations have actually been solved and the boundary conditions. 
inserted. For equation (7a) is not an equation for /; it restricts the lines under 
test by (6), as an equation between x, y, 2, x, y, %, A, in a way analogous to the 
restriction that would be imposed by an equation x = x(t). Only for one particular 
value of h is the world line which minimizes (2) or (6) contained in a set restricted” 
by (7a). Moreover, equations (2) or (6) require the integral to be an extremum 
for given limits in space and time. If we now minimize (9), we obtain a minimum: 
regardless of time limits, as the value of the integral does not depend on time any - 
longer. The minimum of (9) will not in general (i.e. for arbitrary 2) pass the end. 
points in space at the times ¢) and ¢,. Thus, although now the integral in (6): 
has the form fds, its minimum is not that required by Hamilton’s principle. 
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It is thus seen that the analogy between Fermat’s and Hamilton’s principles 
fails, because the lines under comparison in Hamilton’s principle are not lines of 
constant energy which can be fixed arbitrarily, and because the extremum required 
is not an extremum regardless of the duration of transit between local points, 
‘such as would be found by considering the projection of the world line on space 
only. 

As pointed out by L. de Broglie (1930) the problem of generalized optics is in 
fact solved by the principle of least action that has developed from Maupertuis’ 
principle. Following Whittaker (1917) * it may be enunciated as follows: The 


‘equation 
Poti p70 
eos my a ee eee ee tw re 10 
5 [(eagz)ae 0 (10) 


Po, to + Ato 


leads to the equations of Lagrange if the lines under examination are subjected to 
the condition 
L=Xg(OLidgj=hs “harbiteary,, eee (11) 


If now (10) is written 


Pit: At, 
) : (S452) d= Oe “es Wo ere (12) 
0g; 


one can replace x, y, 2, v in the integrand by / and purely spatial derivatives as 
before, allowing at the same time for the condition (11). The line integral is then 
an integral in space only and the open boundary conditions for (12) are satisfied by 
disregarding the ¢ limits. Hence the principle becomes 


Po, tot Aty | 2 


Pi 
OAS. ta Sa ee Awe mee Sere (13) 
Po 
with pe=(lie)eeheLieg,)*. 5 > eee (14) 
With L given by (4), 
p=yimg+(efe)(SA)ys & 3:1) See (15) 


where y is an arbitrary constant, S is the unit vector in the direction of the line, vis 
the velocity of the electron which, as a function of position only, is given by (11), 
Viz. 


—$mv?—eD=h, i.e. v=1/{—(h+e®)(2/m)} ...... (16) 


and A is a constant which can be determined, e.g. by requiring v to have a given 
value at a point of given ®. Then 


p=y[/{—2m(h+e)}+(ele(SA)]. ese (15a) 


If the electric potential is constant, i.e. in the purely magnetic case, v is constant and 
need not be replaced by (16). Alternatively we can express v in terms of (Hp), 
viz. mv =(e/c)Hp, when 


p=yele){Hp} {1+(SA)(Hp)} sane (158) 
which in the purely magnetic case is equivalent to 
pele (SA) p) Oo en eee te (15c) 


* The proof given by Whittaker is especially adapted to mechanics. ‘There are a number of 
‘Jess general formulations of the principle which do not lend themselves readily to the present 
problem. A proof of the equivalence of the principle of least action to Lagrange’s equations, 
‘based on genéral methods of the variational calculus, is given by Bolza (1909). 
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(15 d) is also relativistically correct. Using the notation 0L/dq;=p; where the p; are 
the components of the generalized momentum, the integral in (10) may be written. 


Netdee (Se eta (17) 


Hence it is seen that in general 4=(pS) where p is expressed in terms of the 
coordinates, direction cosines and spatial derivatives by means of (11). 


§3. THE TRANSITION TO WAVE OPTICS: RELATION BETWEEN 
THE RAY AND THE WAVE NORMAL 

Although, as explained in the introduction, the geometrical optics is completely 
described by the differential equations of motion, for the development of the 
physical optics of a system it is necessary to regard the trajectories as rays defined 
by the refractive index and Fermat’s principle. For then the well-known surfaces. 
of Malus associated with families of rays can be established and these in turn are to 
be identified with the wave surfaces of physical optics. This is the transition to 
wave optics. ‘The presence of the term (SA) in the expressions (15) for the 
_ refractive index in electron optics shows that the electron optical properties of 
magnetic fields are those of an optically inhomogeneous and anisotropic medium. 
The surfaces of Malus will therefore not be normal to the rays but they must still. 
be identified with the wave surfaces. 

In an instructive paper, Frank (1933) has shown how the differential equations. 
of a ray in such cases can be expressed by a vector equation (thereby extending the 
method and results of Sommerfeld and Runge (1911)), and has shown how the 
relation between the ray and the wave normal can be determined. The latter 
relation will now be established here directly from Fermat’s integral by extending 
the analysis frequently given for isotropic media (e.g. Forsyth 1927, pp. 229, 
256-262) to the case of anisotropic media. 

In order to link up geometrical optics with physical optics it is necessary to 
consider families of rays, viz. the totality of rays emerging from one point. For 
such a family 

AD, 

Wer \eidt ah eee po TIN ees (18) 
where dr is an element of path along a ray, has a definite value at any point P, in 
space through which a ray coming from Py passes. W expressed as a function of 
this point (i.e. the end point) is called Hamilton’s characteristic function, and the 
surfaces W=constant are called the surfaces of Malus. 

Let us commence by considering a family of non-intersecting lines (in general 
not rays) all emerging from P, and filling the space continuously. ‘They may be 
represented by a family of equations 


w= x;,(u), y=y,(u), EAN Wg Ul RA nis (19) 
where k is a parameter and for all lines x(0) =x) etc. Then each point P, in space 


is defined by a particular set of equations, and a particular value of the auxiliary 
variable u. And, conversely, 


W=| edu, wees (20) 


has a definite value at every point in space, if g is defined by 
Pee! | ways (21), 
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where the refractive index p, originally a function of x, Wz andthe notindependent — 
-direction cosines dx/ds, dy/ds, dz/ds, has become in (20) a function of the 
independent differential coefficients 


x, =dx/du, y, =dy/du, es /d0, ee (22) 
viz. 6 =p(x, V, %, %4/A, y,/A, 2/A), and A is given by 
hNedsdn=s/(a Pty te): See (23) 
Because of its origin the function g satisfies the permanent identity 
= x,(0g/0x,) +.y,(0g/0v,)+2,(0g/02,), =» «sees (24) 


so that dW’ =gdu from (20) becomes 


e a op 
s du+% 50du ha Be Sekt s it Sete 3: (25) 


A Rae! 
dW SAaue aaa: ri ay; 


Hence : 
OW’ /0x =0g/dx,=n,; OW’'/dy=dgl/ey,=n,; OW’ /dz=dg/dz, =n}, 
(20) 
or n’=grad W’, i.e. the vector n’ defined by its components lies perpendicular to 


‘the system of surfaces W’ =constant. 
Now from (21) 


dg = dd + Adp =(p/A)(*ydx, + y1dy, + 2421) 


a inacons wait au c a a # 
+a(Bav+ Ba RH aes Ene *) a(¥) oth ofS 
ant ay t B59 + Be“ \A) © Ayla) VA) * aay \r 


eae (27) 
and 
A(x,/A) = dx,/A — (x4 /A)(xydx, +.y,dy, + 24d2;) ete. 2. (28) 
Thus : | 
7 eh i ema 38) Be aiid Oe ny eons age 
SCE r O(x,/A)\A A(yx/A) r8 O(24/A) r8 
ae Om Dy On Ou Ou 
=pS,— oS, -F (5 105; tS, tage). setegalece, (29) 


In the last expression S,, S,, S, denote the components dx/dS etc. of the unit 
vector along the lines which by (22) and (23) have the values x,/A, y,/A, 2,/A. 

Up to this point the lines x = x(w) etc. have been arbitrary, and in consequence 
also the surfaces W’=constant. Now we specialize: We let the lines be rays so 
that the unit vectors (S,,, S,, S,) along the lines become the unit vectors (R,,, ser 
R.) along the rays and we can identify the auxiliary variable w with the length of the 
* rays measured from the origin so thatA=1. Then n which we write now for n’ is 
given by its components 


= ea a gh 


and W’ ee W now given ne 


Om Ou Ou 
car, + Pvar, +Righ] ese (30) 


-P; 
W= jue Sone (18) 


| 
f 
| 
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the integrals taken along the rays. That means, the surfaces W=constant have 

become the surfaces of Malus. Equations (26) show that the vector n defined by 

(30) is normal to the surfaces W=constant so that (30) gives the general relation 

between the surface normals and the direction of the ray. They are inclined at an 
angle £ given by cos€=(R, nj. By (15), » in the direction of the ray is given by 
: yime 25 (e c(R,A, +R,A, +R.A), 

so that from (30) 

m=>;(moR+t(efc)A),  —= «| 2scce. (31) 

a result given by Glaser (1933). It is seen that n is proportional to the 

‘generalized momentum of the electron.* Hence 


n mzR+(ec)A 


—_ = SESS SO 6 ae eee & (32) 
m \ m*z*+(e c)-A*+2mz(ec)Acosé 
~ where @ is the angle between A and the ray. Then 
ep seine c)Acosh _ / ~ ~(ee)Asind \2 
: V m?z* +(e c)?A?+2mz(ec)Acosé N mv+(e/c)Acosé 
| giving 2 ae He eee (33) 
p (e ec)Asin# 
oe! 
care == c)Acos 5) LS) te ee oem (34) 


Geometrical optics is expressed through the value of the refractive index and 
Fermat's principle alone. The latter assumes that rays can be defined. Wave 
optics is expressed through the refractive index, Huyghens’ principle and the 
wavelength. In general the descriptions are different, in which case that given by 
geometrical optics is wrong. Under certain circumstances, however, they lead to 
the same result and when, but only when, this is true, light can be said to travel 
along the rays of geometrical optics, and Fermat’s principle can be deduced from 
that of Huyghens. Then the surfaces of Malus can be identified with the wave 
surfaces of physical optics and (34) gives the angle between the wave normal and 
the direction of the ray. At the same time the refractive index of geometrical 
optics becomes inversely proportional to the wavelengths. 

This description is based purely on classical physics. To complete the 
description of any electron-optical phenomenon it is necessary to introduce the 
absolute size of the wavelength. Inthe case of anisotropic medium, in the absence 
of a magnetic vector potential, the well-known formula of L. de Broglie 


Io hpnro—chie(iipy Se ee (35) 
is valid. 
As in general MA ee ey es -  AA ees (36) 
equation (154) yields at once 
Aah liat- (SAW Eip) ice Say, . a age (37) 


Here A is the actual wavelength, expressed in terms of the wavelength Ay which the 
electron would have at the same point in space in the absence of A, or at a point in 
space of the same electrostatic potential where A=0. 


* So that curl p always vanishes for a family of rays emerging from a point. Gabor (1945) has 
shown that curl p also vanishes for a family of rays issuing from an extended source, under certain 
restrictions, and hence, that ‘“‘ surfaces of Malus” can be constructed even in this case which has 
no analogy in light optics. "This phenomenon is associated with the existence of the first cross- 


over in front of thermionic cathodes. 
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§4. THE OPTICAL SIGNIFICANCE OF THE ARBITRARINESS 
OF THE MAGNETIC VECTOR POTENTIAL 


The detailed discussion of the origin of given above shows clearly the extent 
to which the mapping of ,. throughout space is arbitrary. The expressions (15) © 
as written down are unique. An arbitrary additive constant in L would have 
dropped out. Therefore the only freedom left in assigning a numerical value to 
yu arises from an arbitrary additive constant F included in ® and the gradient of © 
arbitrary scalar % included in A: A=A,)+grad;. Hence in general 


w=y{[—2m(h + e®, + F)]} + (e/e)(SA) +(e/c)(Sgrady#)}, ...... (15d) 


where the first term denotes the isotropic contribution to mw arising from the 
electrostatic potential and the second and third’‘terms the anisotropic contribu- 
tions. # combines with the arbitrary constant of integration and produces thus 
no new effects at all when / is determined by the boundary conditions. ‘There is 
no general term k/v (k an arbitrary constant) as proposed by Opatowski in 1943. 
The term containing grad % to which he has rightly drawn attention is, however, of 
considerable interest. 

Electrodynamic considerations require only that curl A reproduces a given 
field H everywhere. But (SA) occurs as an additive term in the refractive index 
and it was seen in §2 that any expression for the refractive index is subject to the 
condition that it is continuous and free from singularities. Hence (SA) and there- 
fore A must satisfy the same condition if it is to play its role in electron optics. 
These conditions are just those for the validity of Stokes’ theorem (Kellogg 1929, 
p. 91). The further condition of the continuity of the first derivative need not be 
specially considered in physics as it can always be satisfied by a physically 
insignificant alteration in the vector potential. : 

Hence it appears that the vector potential in electron optics must satisfy 
Stokes’ theorem, and this is the only valid restriction which must be imposed on it. 
In consequence A cannot in general be chosen so as to vanish with the magnetic 
field. An example will be given in the last section. 

In the electromagnetic theory it is customary so to choose A that for static 
fields divA=0. 

This condition is frequently applied to the vector potential in electron optics, 
as, for example, in the book by Zworykin et al. (1945). It is readily seen, however, 
from the development given, that in this connection such a condition is entirely 
arbitrary. Further, it is not useful as it does not fix the value of A uniquely, for the 
value of y can still be chosen to be any solution of Laplace’s equation div grad ¢ =0, 
and is therefore not uniquely defined. The remaining arbitrariness in A, or 
grad ¢ if an Ap satisfying Stokes’ theorem is given, is still considerable. 

As the arbitrariness of the vector potential does not enter Lorentz’ equation it 
cannot produce any observable effects in geometrical optics. This can, of course, 
also be shown by a consideration of Fermat’s principle with the value for p given. 
But equation (34) shows that the inclination € of the wave surfaces to the rays 
depends on the absolute value A of the vector potential and on its inclination 6 to 
the rays. As the addition of a term grad#% to Ay leaves the rays unchanged but 
changes A with respect to both direction and abel value, both A and @ depend 
on the arbitrary ys. The form of (34) shows that in general the change of A is not 
compensated by the change of 6, so that in general the addition of a term grad % to 
A) results in a change of &, and thus, as the rays remain unaltered, in a change of the 
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wave surfaces. Hence by a proper choice of 4 an anisotropic medium can be 
made isotropic locally, and an isotropic medium can be made anisotropic, in an 


infinite variety of ways. 


It is seen, therefore, that the extra term (e/ce)(S grad) profoundly alters the 
characteristics of the electron waves, but it will now be shown that it does so insuch 
a way that there are no consequences of this change that could be observed in any 


_ experiment. The basic problem of physical optics is the evaluation of the optical 
_ path lengths of pairs of interfering rays connecting given points. Apart from 


questions of intensity, the physical phenomena observed depend only on the optical 
path difference for such pairs. Now the optical path length along a given line 


- connecting a point 0 and 1 is given by 


J c ds= | [mo+(SA)]ds, see (38) 


for a particular choice of A, viz. Ay. For another value of A given by A, + grad ys 
1 

the optical path lengths will be different from (38) by | (S grad i)ds =, —p. Ass 
; 0 


is single valued this difference is the same for each ray connecting the two points 
considered. Thus the physical result of any electron-optical calculation and in 
particular any interference effects is independent of how A is chosen but the 
detailed picture of the relation between rays and wave surfaces is not. The 
particular set of wave surfaces belonging to a family of rays and the inclination of 
the wave surfaces to the ray at that point will vary according to the choice of A. 
This is compensated, however, by the change in the wavelength given by equation 
(37). By inserting reasonable numerical values for A and (Hp) it is readily seen 
that the factor of A) may be considerably different from unity, so that the wavelength 


- of anelectron beamis largely arbitrary. 


Since the physical and geometrical character of a solution of an electron- 


optical problem is independent of the particular choice of A, one may expect to 


find a solution of .any electron-optical problem in terms of uniquely defined 
quantities. Such a quantity is the magnetic field H, and it is easy to see how 
problems can be solved in terms of H rather than A, for the solution requires only 


_ that it should be possible to evaluate the difference d in the values of J uds along: 


any two paths between the object point and the chosen point of the image space. 


- With p given by (15) 


a=(ele)y| i ' (Hp)ds — | “(Hp)ds' + | (SA)ds — if (SA), sek (39) 


_ where ds denotes elements of the first and ds’ those of the second path. Now,, 


| “(SA)ds— { (S’A)ds! = $(Ads) = |e eee (40) 


where the last integral but one is taken round the closed path formed by the two 
lines under consideration and the last one arises from it by Stokes’ theorem and is 
the integral of the normal component of H over any surface bounded by the two. 
lines. Hence 


d=(elc)y (| (Hide | pds + | [H.ae), er (41) 
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In general (Hp) is a function of the coordinates; in the absence of an electric 
field (Hp) =const., when 


; 
d= (ele ({ zy -= | Hd) eee 42 
(e/e)r(He) (4 ca || (42) 
and the order of interference for the two lines comes to 
1 1 
(oh + zaps | [Hado) er: 43 
ale A+ cai | Hato fe 


where /, and J, are the length of the two lines joining Py and Py. 

An elementary illustration of (43) is the deflection of a parallel bundle of 
electrons by a thin magnetic field of strength H and thickness b. If the bundle is 
deflected by a small angle ¢, two rays d cm. apart acquire a path difference ¢d which 
by (43) must be compensated by {1/(Hp)}Hbd = Hbde/mcv as the deflected beam ~ 
is still of zero order of interference. Hence d= Hbe/mcv which is a well-known 
formula. 

§5. SPECIAL PROBLEMS 

It remains now to show how the refractive index given leads in fact to correct 
solutions of electronic problems, and to a consistent description of electron- 
optical phenomena. ‘ 

A. Differential equations for trajectories: For che solution of particular — 
problems it is convenient to mark out one of the three coordinates q,, qo, d3, e.g. x if — 
the gq; denote cylindrical coordinates x, 7, 6. ‘Then generally ds=Adq,, e.g. for 
the case mentioned ds = 1/{1 + (dr/dx)? +.7°(d0/dx)*}dx. ‘Then Fermat’s principle 


G1142'ds* 
5 pid. = 0 Ae le eee (44) 
1 We Ws 
leads at once to two independent differential equations for trajectories 
4 AwA) _AwA)_9, @ _AwA)__ Awd) _ a 
dq, O(dq2/4q;) 0g. ; dq, 0(dq3/dq;) 09s eer ; 


The expression for A is greatly simplified if the field has symmetry, such as axial 
symmetry. For this particular case «A = mvA + (e/c)Agr(d0/dx) and the equations 
(45) lead through a short calculation to the fundamental equations of electron lens 
optics (see, for example, Zworykin et al. 1945, p. 503). 

B. It may be of interest to illustrate the discussion of the wave surfaces in 
electron optics given by reference to a simple two-dimensional example. 

Korsunsky (1945) pointed out that if a parallel bundle of rays in a plane were 
incident on a uniform magnetic field perpendicular to the plane, having a circular 
boundary so that H = Hy for 0 <r <a and H=0 for r>a, then these rays will be 
focused. ‘The field acts asa focusing prism. If the “‘ momentum” cmv/e =(Hp) 
is such that (Hp)= Hoa the prism gives a perfect focus over its entire aperture; 
this focusing action is illustrated in Figure 1. The rays are straight lines outside 
the field and circles of radius a inside it. In this case the optical paths of all the 
rays arriving at the point I should be equal, and that they are is easily tested using 
equation (42). 

Any ray A,BI of radius (in the magnetic field) a, can be distinguished by the 
parametric angle €, and any point P on the ray by either the angle 0 or the angle ¢ 
constructed as shown in the figure. @ increases outside the field from 0 to € 
corresponding to the point B where the electron enters the magnetic field, and 
inside the field from € to 7/2 at I, while ¢ increases inside the field from zero to 
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2$,=7/2—€ at I. The angle ¢, corresponds to the closest approach of P to the 
axis O of the magnetic field. 

Then, using equation (42) but leaving out the constant factor (e/¢)y(Hp), the 
| optical paths A, BPI and Aol coming from infinity differ by the difference of their 
| actual lengths plus Hoo/(Hp) =a/a where o is the area between the arc BPI and the 
edge of the circle. The actual lengths are /,,+a(1—cos£)+2ad, and /,,+4a 
respectively, where /,, is the length of any ray from infinity to the line MM’. 
| The area o is 


2(3a°24, —a* sin $, cos $,) = a*(2, — 2 sin ¢, cos d,) = a(2ad, —acos €). 
| ‘Thus it is verified that the optical paths are equal. 
In order to construct the wave fronts a vector potential describing the given 


' field is required. It is readily seen from the expression for curl A in cylindrical 
_ coordinates that 


A,=A,=0; Ay=37rHy, Nxrsal A,=@H/r), a<r wt... (46) 
H = 
I 
a. Figure 2. 


Figure 1. Figure 3. 


satisfies this condition. If A, had the value given in the inner region but vanished 
in the outer region, an infinite field of opposite sign at the boundary would be 
indicated, owing to the discontinuity of Ay. Hence the need for a continuous Ay. 
‘The field discussed here is the limiting case of a field continuously falling off 
from H, to 0 in the neighbourhood of the boundary. It is easily seen that the 
simplification introduced by letting the field fall. off discontinuously has no 
significam effect on the results given here. Consider now the optical path length p 
of the ray A,BP coming from infinity to a point P outside the field, and let P be 
given by its distance « from MM’ and the angle 0=0". 


p= (* eS {4 5) eae al.) 


d= 


aH, 


sin 6 ds, 


since A is always perpendicular to 7 the radius vector from 0. With sin 6ds=rd0 
the optical path considered becomes 
Dia oat Gis ety VO Ae ee (47) 


2-2 
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If P is inside the field, the part outside the field is given by (47) with «=a(1—cosé) — 

and 6’=é. For the part BP in the field (AS) =47H)sin(6+¢4); ds=ad¢, so that _ 

the optical length of the arc BP is given by af+ 3 [ rsin(6+¢)d¢. Since the © . 

chord BP of length 2asin¢/2 makes an angle ¢/2 with the direction A,B, 
rsin@ =asin€é+2asin?4/2; rcos@=acosé—2asin¢/2 . cos ¢/2, 

then rsin(6+¢)=asin(é+¢)—2asin?¢/2. Thus 


4 [' rsin (6+4)d6 =} —acos (€+¢)+acosé+asind—ad}, 


and the total optical path from infinity to the point P is 
p=l,,tattafé—cos€é+sind+d—cos(E+¢)}.  -- eee (48) . 


Notice that if P goes to I, i.e. ¢ goes to 7/2 — & this expression becomes /,,+a(1+7/4) 
a constant value indépendent of €, confirming that all the rays Ap, A, etc. have the 
same optical path to I. 

Without restriction in generality, we may write /,=0 and a=1, when the 
wave surfaces (lines of constant p) outside the field can immediately be constructed 


from (47) in terms of x and 6’. Inside the field, p(¢) may be plotted from (48) for — 


a few rays given by their parametric angle €, when points of given p are constructed ~ 
in terms of € and ¢. . 

The wave fronts are shown in Figure 2. There is no discontinuity at the © 
fieid boundary and the medium is strongly anisotropic outside the field. 
According to §4 the wavefronts shown are not unique. 


One may ask if the anisotropy outside the field could not be avoided by an 


alternative value for A which also reproduces the field given, such as 
A,=A=0,  Ay=t7h, —3¢H, 7, Osraa@: A,=02 aan 


So chosen A vanishes with H outside the boundary of the field, yielding here an 
isotropic refractive index, and a wavelength given by de Broglie’s formula. While 
this A correctly describes the given field, it does not satisfy Stokes’ theorem since 
A, becomes infinite at the origin. | 

In fact, the vector potential (49) can be derived from (46) by the addition to the 
latter everywhere of the vector whose polar components are A,=A, =0, 
A, = —3a@°H,/r, that is, of the vector grads where = —}a?H,6. The new term 


nL 
would add to the integral Ee weds considéred above an amount—4a?H,(6,,—6;). 


This has one value for all rays that pass below the origin O and another for all 
rays passing above it. If @ is measured as indicated in Figure 1 and €, denotes the 
ray passing through the origin, the effect of the new term is to add the constant 


“I 
value aH 7/4 to the integral | weds for all rays €,>&>7/2 while it adds 


: =e} 

— @H)(37/4) for all rays €,>€>-—7/2. This means that the oncoming electrons: 
cannot all be described bya single plane wave but that there are two plane waves 
necessary that cannot be related in phase. This effect is associated with the 
breakdown of Stokes’ theorem for grad 0, the 6 component of which is infinite at 
the origin. 

It is readily seen that no vector potential which satisfies Stokes’ theorem will 
remove the anisotropy of the whole space outside the field, so that a consistent 


Refractive Index in Electron Optics 21 


_ optical description appears to be impossible while the whole field-free space is 


anisotropic. The irremovable anisotropy of the field-free region as a whole 


_ emphasizes the fact that the electron-optical refractive index contains the vector 
| potential and not the magnetic field strength. The formulation for the optical 
¢ path difference given in equation (42) which does not contain A does not alter this 


situation, as a surface integral over H is equivalent to a vector potential. One 
might therefore expect wave-optical phenomena to arise which are due to the 
presence of a magnetic field but not due to the magnetic field itself, ie. which arise 
whilst the rays are in field-free regions only. Consider now an arrangement as in 
Figure3. Odenotesa point source of electrons which is focused by the lens M at 
the point P. Through the pair of slits separated by / a set of interference fringes 
will arise so that the distance of the mth maxitnum from P is given by d=bAn/l. 
If, now, a magnetic flux is established normal to the plane of the paper through 
the area a, then, according to (43), the order. of. interference at any point of the 
focal plane is changed by 


ere A 
iF \o(Hp) 


Thus, a flux of 3-9 x 10-7 gauss cm? is required to change the order of interference 
by 1, and half of this flux will change the maximum at P to a minimum. 

It is very curious that equation (50) associates a phenomenon observable at 
least in principle with a flux; one expects a change of flux, but not steady flux, to 
have observable effects. ‘The effect has, however, a certain analogy in the existence 
of a permanent current in a superconducting ring due to a magnetic flux through it. 


N | | Hideo, or with (35) by N-=(e/ch) | | Hideo. Br, (50) 
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The Disturbance near the Focus of Waves of Radially 
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ABSTRACT. Lommel’s problem, the disturbance near the focus of waves of uniform 
amplitude, is generalized to consider waves of radially non-uniform amplitude. Lommel’s 
numerical results are extended, and his U and V functions, tabulated on a new basis, 
together with some new related functions, are employed to provide extensive intensity 
distribution curves for the diffracting patterns associated with waves of both uniform and 
non-uniform amplitude. 


§1. INTRODUCTION 


was derived by Airy (1834). It has the form 
1 S425, (2) 24, sige i 


where J is the intensity at points distance p’ from the centre of the diffraction 
pattern, and 


dee distribution of intensity in the geometrical focal plane of a spherical wave — 


2=2nrp /AR. Ape ees) 


The radius of the circular aperture which limits the wave is 7, and R is the distance 
between the limiting aperture and the focal plane. . is the wavelength of the 
disturbance in the medium in question. - 

The expression (1) accounts for the form of the diffraction pattern, consisting 
of a bright central maximum surrounded by alternate dark and bright rings, which 
is observed as the image of a distant star in a well-corrected telescope. The 
dependence of the scale of the diffraction pattern on the relative aperture and the 
wavelength is shown by (2). 

In deriving (1) it is assumed that the angular semi-aperture is small, and that 
the wave is of uniform amplitude. ‘The former restriction has been removed in an 
earlier communication (Hopkins 1943), the error in the approximate formula being 
shown not to exceed a few per cent for « <30°, where « is the angular semi-aperture. 

With the approximation implied in Airy’s formula, Lommel (1886) considered 
the intensity distributions in planes other than the geometrical focal plane. 
Development in series of the integrals concerned led to the introduction of the 
Lommel U and V functions. 

Itis proposed here to remove the restriction of uniform amplitude from Lommel’s 
problem, and to investigate formulae expressing the disturbance near the focus of 
a spherical wave showing a radial non-uniformity of amplitude. It is found 
necessary to introduce both the Lommel functions and new related functions. 


These functions have been tabulated and used to provide extensive numerical 
results. 


— 
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82. LOMMEL’S PROBLEM FOR WAVES OF RADIALLY NON-UNIFORM 
AMPLITUDE 


In Figure 1, O(X, Y, Z) is a system of rectangular coordinates, O being the 
_ centre of an approaching spherical wave. P(X, Y, Z) is a point on this latter, and 
) YOZ is the geometrical focal plane. The wave is limited by a circular aperture 
| of radius 7, and OP=R is the radius of curvature of the wave. 
For small angular apertures one can write pdpd¢ for the area of an element of 
| wave:at P, where (p, ) are polar coordinates of P, having origin on OX and initial 
line in the plane ZOX. The disturbance at P’(AX, Y’, Z’), a point in a plane of 
focus a small distance AX from YOZ, can be written 


ge [ ; S BOAT dijadteh Oe (3) 


where d= PP’, Ais the amplitude at P, Tis the period of the disturbance, and ¢ the 
time variable. 


The distance d is given by d?=(X —AX)?+(Y—Y’)?+(Z—Z’)? or, AX, Y’ 
and Z’ being small, d?= R?—2pp' cos(¢—¢’) —2X AX, where (Y,Z), (Y’,Z’) 


xX : : 
H 
Ss DERE CENTRE oO! 
PEG) Pp WAVE CENTRF © 
z 
Se 
° 
' o+ 
ew PYCAXsY.2) 
Figure 1. Figure 2. 


have been changed to polar coordinates (p, ¢), (p’, 4’), the latter referring to the 
point P’ in the plane of focus through O’. The coordinate X is given by 

=(R? —p?)? = R—(p?/2R) with the required accuracy. ‘Thus 
d=(R—AX) —(pp'/R) cos (b—¢’) + (p?/2R?)AX. vexens 4) 


Substitution of (4) in (3), and the expansion of the argument of the sine in the 


integrand in terms of those parts dependent upon, and those independent of, the 
variables p, ¢ gives 


= 1/(C? +S?) sin 2n(t/T —(R-AX)/A+e), 
where tan e = C/S and 


Z ; aoe 
Sa ff ASS AB cos(b—$')— Tee OX} papas, 
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Thus the intensity at P’ can be expressed as the squared modulus of the integral - 

Cus= i) fr Aexpi{(2mpp’ /AR) cos (6 —¢$") —(mp?/AR*)AX} pdpdd. ...... (5) 
If now the amplitude A varies radially, one may write (5) 


C+iS= |. Aexp{ = s(mpi/AR) AX} pdp ir exp {1(2pp'/AR) cos(¢—¢')}db — 


ar IE "Aexp(—iyl2\(e/r)*Joa)ede, nee (6) 


where y/2 =(zr?2/AR?2)AX, and x=(27/A)(pp’/R). The significance of these will 
be seen from the following. In Figure 2, HP isa section through the spherical 
wave-front, and HQ a section through a sphere centred on O’. The marginal 
difference PQ measured in wavelengths, is 


W =(PQ)/A=r?AX/2AR?, of .y/2 = 2a We 7 Oe atone (7) 
It denotes the amount out-of-focus, independent of the aperture. 
Again, let z be the value of x when p=, then p=rx/z; pdp=(r?/z*)xdx; and 
z=2mrp'/AR, which is the variable defined in (2). It denotes the scale of the 


diffraction pattern independent of the aperture and wavelength. 
The integral (6) can now be written 


2 
Cs 
a 


~i(y/2)(2]/2)?]So(x)xdx, wea. (8) 


which is Lommel’s integral if A is taken to be constant. Write for the amplitude 
A=Q+P(p/r)?. Then A will have the value Q at the axial point of the wave-front, 
and the value F=Q + Pat the edge of the aperturep=r. ‘The amplitude decreases 
or increases on going from the axial point of the wave-front according as P is 
negative or positive. Equation (8) now becomes 


Casa =. {0+ P(: y hexp | = (=) * [roleva, 
“r(o +2iP — mf | exp | - iz (=) "| wel (oc dae ee iy 


The integral to be evaluated in now simply that of Lommel’s original problem.* 
This is expressed in two forms, suited respectively to the regions s<y, s>y. 
These‘ results are (Gray and Mathews 1922) 


Cae aN Reged 
2 [‘exp| -3 (5 =) | stola)ax= 5 (sins + Vsin by ~ V,c0s 4y) 


2 2 5 aa 
ah Gone —V,cos4y—V,sin y) Sie aS Cera (10) 


or 


Cie 


for the region z<y, and 


Dee 
z vexp| -13 (2) eee sindy) 


+42 (U, sin by=Ucosky\ a eee (11) 


* This step shortens the analysis. I am indebted to a referee of the paper for pointing it out. 
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_ for the region z>y, where 


leg) 


U, ae a 1)?(y/z)? 72°F n.05(8)3 ka = py ( “#7 1)?(a/y)"t?P Sn 4on(8)- er Speen (12) 
= p=0 


The results (10) and (11) immediately yield expressions for (9) suited to these 
two regions. ‘Thus, omitting the constant factor zr”, we have, for the region z<y, 


; Pee 
C+iS= (o+21P 2 =) (sin + Vysinsy— V, 008 49) 


NO NES 
+1(Q+2P 5)" (cos — V,cos4y—V, sin iy), 


or, with some arrangement of terms, 


C= FT + P(2/y)(l',—%); S= hie PZ/y\U kin) ieee (13) 
where 
r vay Ba hele ‘1 } 
oi SA +Vosindy— V; cosy), 
9 snk EN ee a Sat ath a Pa MS (14) 
a= 5 (cos 5 —V,costy— V; sin dy) : 
and 
a ees? 
r= 5{(F - 5) sin = + Yo cossy+ Y, sindy}, 
a Depa ? of Se ag (15) 
Bie \\ Oyo 5) cos 5 + osinsy— cosy}. 
‘The Y functions and the X functions are defined by 
il 0U,, fee) , y n+2p 
X= 5 Se = E1900 +20) (2) Ty sal | 


: Soy Sean (16) 


o a\nt2 
Yum — ype = 3 (= 1)(0429)(S) Tsar) | 


For the region z>y 


C+iS= (0+2P 5) <(U, aa U, sin $y) 


+i(Q+2iP) (U, sinty — U, cos $y); 


expressing C and S in the form (13), one must now write 
[T=(2/y)(U, costy+U,sinky); X=(2/y)(U,sinty—U,coszy) ...... (17) 
and 
=(2/y)(X, sin }y —X,.cos $y); U4 =(2/y)(—X1 cos 2y —Xsin 2y). 


It will be seen that in the case F = QO (i.e. P=0) the above formulae degenerate 
to Lommel’s results for waves of uniform amplitude. Further, the expressions 
M13); together with (14), (15), or (17), (18), provide a solution of the present problem, 
for the intensity is given by J(y, z) =C?+.S?, and can be calculated at once given 
tables of the U, V, X and Y functions. 
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§3. GENERAL FEATURES OF THE INTENSITY DISTRIBUTIONS 


As one would expect, the intensity is independent of ¢’. If be replaced by 
—y, the expression /=C?+ S? remains unchanged, although C changes sign. 
Thus the intensity distributions, considered as functions of 2, are identical in 
planes equidistant from the plane y=0. Hence, as in the case of waves of uniform. 
amplitude, the spatial distribution of intensity is symmetrical both about the 
geometrical focal plane and about the central line z=0. 


(i) The Geometrical Focal Plane 
Putting y=0, one obtains for the geometrical focal plane 
C= F(2J,(2)/2]— P[4J,(z)/z2],  S=0. 
The intensity is thus 
I(0, 2) = {F[2J,(2)/2] — P[4Jo(2)/27]}?, 


which yields the Airy formula (1) for F=Q=1, P=0. The intensity at the © 


geometrical focus is 1(0,0)=(Q+4P)? in the general case. 


(ii) The Central Line z=0 
Taking the limits of (14), (15) as z-0, one has 
C= F(2|y) sin dy + P(2/y)*(cos $y ~ 1), 
$= F(2y)(1—cos 4y) + P(2/y)%sin dy) —2Py. | 
A zero of intensity requires that C, S be simultaneously zero. ‘That is, since: 
F=Q+P, 
Fsinty+(2P/y) costy=2P/y; —Fcosty+(2P/y)sinsy= —Q, 


sin ‘i 


or [P+ QP) cost +9= 6) 


Squaring, and adding these equations, F=Q. ‘Thus there are no zeros of intensity 


along the central line, except in the case of waves of uniform amplitude. 
Writing F=Q-+ P in (19), we find for the intensity along the central line 


sin + 2 inly 
I(y,0)=(Q+4P) (=) +P? (<) (= sos ty) ee ON. (20) 
ay y. ay 

which is easily interpreted since }y=7W as in (7). 

If Q=F=1, P=0, (20) gives I(y, 0) =(sin}y/4y)? as the intensity distribution 
for waves of uniform amplitude. ‘The zeros occur at integral values of W, 
excluding W=0; the maxima occur at the roots of tanm7W=7W. The maxima 
of (20) again appear at the roots of tana W=7W, and the minima at the roots of 


20(0 + P)\(rW)? tana W — P°(nW—tanaW)=0. 


The minima are, in fact, displaced towards the positions of the maxima—degener- 
ating, as the above equation shows, to points of inflection for the cases Q =0 and 
F=Q+ P=0 respectively. 

At points along the central line corresponding to integral values of W, the 
intensity is zero for waves of uniform amplitude (P=0). For non-uniform waves 
the intensity at these points is [(47W, 0) =(P/27W)?. 


oo ec earmar 
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Thus, at the first “dark point”’ the intensities, in different cases, are as follows : 
P= 1-0: 31(47,0) =0-025, £(0,0) = 0-250, 
P=—0-4:  I(47,0)=0-004, (0, 0) =0-640. 

The first of these cases is that of a wave whose amplitude decreases from 1-0 to. 
zero on going from the centre to the edge of the wave. The intensity at y =47 is 
10% of that at the geometrical focus. The second case is that of a microscope 
objective, for which the non-uniformity of amplitude has been calculated using 
Fresnel’s formulae in a manner shown earlier (Hopkins 1944). At the first 
‘dark point” the intensity is 0-6°% of that at the geometrical focus. 

Consider now two waves, in one of which the central and marginal amplitudes 
are O=a, F=6; and in the other Q=f$, F=«. In each case (20) becomes 

_ [a+P\? /sincW\2 c—B\2 (1. \2/sinnW ) 
ez) Go) ar eh 


7 


the interesting result that the intensity distributions along the central line are the 
same. However, considered as functions of z, the two cases differ widely. 


$4. TABULATION OF THE U, V AND X FUNCTIONS 
Lommel has provided values of the functions U,, U,, Vo, V, for y=z7, 


Pa,.... 10m, and 2=0.1,.... 12. For problems in the diffraction: theory of 
spherical aberration, using a method developed earlier (Hopkins 1945), it is 
desirable to cover the ranges y= — 0 to + ooandz=0,1,....36. The following 


method of tabulation was therefore adopted in each case. Writing A =y/z, and 
B=32/y, the definitions (12) and (16) give 
£0 n+2p ee N+2p 
y= 2(-0(S)  Inamleds y= Z(—IP(B) 4+ 20 nsnl®) 
pias. sone >, U-= VV X=. 

For even values of n, the functions do not change sign with the argument A or 
B, while the sign is reversed for odd values of n. A table showing U,, X,, for 
ai 01)-05,95,,. 1-0 and 2—0, |,.... 36 willthus-eover U,, V,,; X,, and Y,, for 
the required ranges. 

The tables were computed by machine, using five-figure values taken from the 
tables of Bessel functions given by Gray and Mathews (1922). 

The functions reduce to simple expressions when A =0 or d=1. ‘These were 
used to check the terminal values, intermediate ones being checked by plotting the 
first differences as functions of A. ‘Three places of decimals were found to be quite 
adequate for the present problem. 


§5. INTENSITY DISTRIBUTION CURVES 


In the case of a light wave originating from an axial point source, the energy 
transmitted along marginal ray-paths through a lens system may differ greatly 
from that along the axis, since here the incidence is always normal, whereas along 
marginal rays large angles of incidence may be encountered. It is not difficult to 
show that the variation of amplitude with aperture should be a linear function of 
the square of the aperture. This is the case treated above. 

In an F/2 photographic objective, for example, the amplitude at the margin 
was but 70° of that at the centre, while in the case of the microscope objective 
referred to above it was but 60%. Thus it seems reasonable to assume that the 
case Q=1-0, F=0-5 would represent the very extreme likely to be met with in. 


28 H. H. Hopkins 


practice. On the other hand, an artificial “inhomogenizing”’ of amplitude is q 
limited at the cases Q=0-0, F=1-0 and Q=1-0, F=0-0. These cases have | 
therefore been computed in order to investigate the importance of the effect in ) 
-common practice, and its theoretical limits. ‘The curves are all shown on a scale 
such that the intensity is unity at the geometrical focus. 


(i) Uniform Amplitude. Q=F=1 
In each figure the intensity curve for this case is shown as a broken line. In 
-each plane of focus the intensity is computed for z=0, 1, .... 24, thus extending 
Lommel’s results, which are given for values of z up to 12. 
It will be seen that the out-of-focus distributions approach the “‘ Airy disc”’ 
(shown as the broken line) for large values of z. It is not difficult to show this to be 
required by the above formulae. 


(ii) The Central Line z=0 


The intensity curves for three cases are shown in figure 3. The marked 


INTENSITY 


AMOUNT ‘“OurT - oF ~- FOCUS” W = WAVELENGTHS. 
SS NESSarasnds soemeiewcsssos 


Figure 3. Intensity along the line z=0. 
—— —— Q=F=1-0. 
———>—  Q=1-0, F=0-°5 (or O=0:5, F=1-0). 
——— 9=0°0, F=1:0 (or O=1-0, F=0-0). 


features that appear as the non-uniformity of amplitude increase are the dis- 
appearance of the ‘‘ dark points”’ and the increase in depth of focus. This latter 
is further exemplified in the curves of (iv). 


(i) Non-uniform Amplitude. Q=1-0, F=0-5 


The intensity curves in this case, those of a wave whose amplitude decreases toa 
value at the edge equal to half that at the centre, are shown in Figures 4 to 12. It 
will be seen that the differences between the forms of the light distributions curves 
of this case and those of a wave of uniform amplitude are not very significant from 
a practical point of view. It is therefore justifiable, in the case of almost every 

-actual lens system, to assume the emergent wave-fronts to be of uniform amplitude 
so far as the form of the light distribution curves are concerned. 
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Figure 4. Figure 5. Figure 6. 


- Plane of focus: W=0. Plane of focus: W=H. . Plane of focus: W= re 


Figure & 2 Figure 9. 


Plane of focus: W=1. ; Plane/of focus: W=14. 


\ 


9 4 8 12 16 20 24 ° 4 8 12 16 20 24- 


Zs j ; Zz 
Figure 10. Figure 11. Figure 12. 
Plane of focus: W=1. Plane of focus: W=12. Plane of focus: W=2. 
O- Ht) —— Q=1, F=4 (on 1-7775 x scale). 
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Figure 13. Plane of focus: W=0. 
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Figure 14, Plane of focus: W=}. 
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Figure 15, Plane of focus: W=#. 
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Figure-16. Plane of focus: W=3. 


O2Fat Oa rads 
QO=1, F=0 (4xscale). —— Q=0, F=1 (4xscale). 
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Figure 17. Plane of focus: W=1. 
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Figure 18. Plane of focus: W=1}. 
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Figure 19. Plane of focus: W=1}. 


Figure 21. Plane of focus: W=2. 
——— Q=F=1. — Q=F=1. 


OQ=1, F=0 (4x scale). —— Q=0, F=1 (4xscale). 
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(iv) Non-uniform amplitude. Q=1-0, F=0-0; Q=0.0, F=1-0 


These are the extreme cases, those of waves of amplitude decreasing to zero at 


the edge and centre respectively. The curves are shown in Figures 13to21. In 


the case Q = 1-0, F=0-0, the tendency is for the central maximum of the diffraction 
pattern to be broader than in the case of uniform amplitude, accompanied by a 
virtual suppression of the ring, system. In the case Q=0-0, F=1-0, the central — 
maximum is narrower and the ring system brighter and showing greater contrast. 

It is thus seen that artificially produced non-uniformity of amplitude can 
significantly modify the light distribution curves, although for most ordinary lens. 
systems it is quite justifiable to assume uniform amplitude in evaluating the 
diffraction integral. 
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ABSTRACT. ‘The results of a spectroscopic study of the white or coloured flames from 
typical pyrotechnic compositions are reported. Most of the visible light from these sources 
comes from a limited number of atomic or diatomic emitters. The effects of secondary - 
emitters on the dominant wavelength, colorimetric purity and relative intensity are con- 
sidered in general terms and some suggestions for improving the quality of the flames 
are put forward. 


| §1. INTRODUCTION 


YROTECHNIC “compositions” are used in certain circumstances, for 
Pp instance in signalling, to produce coloured or white light of high intensity 

and moderate duration (of the order of seconds). They consist of compressed 
mixtures of finely divided substances which, when ignited, react vigorously to 
produce luminous flames. Such systems are clearly very complex, and we shall 
consider only two main aspects: first, the spectroscopic characteristics of the 
flames, and second, the relation between these properties and the visual sensation 
that is produced by the flame. 

The considerable energy required for the production of light from such a flame 
is alrnost always derived from an oxidation reaction. For such flames, the com- 
position contains (i) a source of oxygen, e.g. a chlorate or nitrate; (ii)a fuel, such asa 
carbohydrate or a metal like aluminium or magnesium; (iii) a substance which will 
impart colour to the flame, e.g. cuprous chloride; and possibly (iv) small amounts of 
substances to protect an electro-positive metal. to produce gas to increase the size of 
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the flame, to act as a binder for the compressed composition, or to modify the rate of 
burning. Sometimes, of course, it is possible to choose a substance to fulfil 
simultaneously more than one of these various functions. These general remarks 
are illustrated by some specimen compositions given in the next section. 


§2. SPECTROSCOPIC PROPERTIES 


The spectra of the flames were photographed in a first order of a 2°4m. concave 
grating in an Eagle mounting, which gives a nearly linear dispersion of about 
74a/mm. Ilford ‘Special Long Range Panchromatic”’ plates were used. The 
spectra were normally examined over the range 3000 to 8000. 

Altogether a large number of different compositions were studied: those 
listed below (taken from Weingart 1943 and Davis 1943) are typical examples. 


Blue es Green 5 
Potassium chlorate 5 Se OS Bariuminitrate,.  . = ..50 
@opper carbonate - - . 15 Barium chlorate he eee otal) 
Mercurous chloride Mi-ptcw wD) See One| acm come -aeerr in a ee) 
Shellach. wasn. > ee 12 

Yellow Red 

Potassium chlorate ee 55 Potassium chlorate ere 
Sodium oxalate MMR eee) Strontium nitrate . . . 65 
Shellacwe. =e. 6 208e, SA $5 Shell acura ea enema () 
White White 

Niigioatianiuten 6 Wg) 6 apm Via OTe SULA en eee) 
Barun mithate ss soa) e875 Barium nitrate . . . . 45 
Sulplaptbes eet. Ra ates wen 3 Paraffin wax Pen) Tita te its) 


Most of the visible light from these sources comes from various atomic or 
diatomic emitters. It will be best to consider the different types of composition in 
turn. 


Blue flames. | 

These rely for their colour on emission of the spectrum of CuCl. Most of the 
light comes from the D and E systems in the range 4200-4600 a. (cf. Pearse and 
Gaydon 1941); the blue-green systems B and C, and the green system A, are much 
weaker (even allowing for the drop in sensitivity of the plates in this region). 
Potassium chlorate and perchlorate are often used both as oxidizing agents and as 
sources of chlorine for the formation of CuCl in the flame; this leads to strong 
emission from the K resonance lines at 7665 and 7699 a. ‘The blue composition 
given above contains mercurous chloride, which, however, can have no function 
except to provide chlorine (Wieland 1945); equally satisfactory compositions may 
be prepared which contain none of this substance. 


Green flames. 

BaCl, which gives rise to a well-known band system at 5050-5350 a., is the 
obvious choice for a green emitter, and it is in fact generally used. ‘The spectro- 
grams of such flames usually show, in addition to this system, band-systems of 
CaCl and of SrCl in the orange and red regions of the spectrum, which evidently 
arise from small amounts of impurities in the original barium compounds, and they 
invariably show the extensive band-system of BaO (4000-8000 .), There is also 
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strong emission from a near infra-red system of BaCl (Barrow and Crawford 1946) | 
which, however, lies beyond the visible region of the spectrum. Other discrete 
bands in the region 7000-7600 a. have not yet been identified. 


Yellow flames. | 

When sodium is introduced into a flame such as that of hydrogen burning in air, 
the emission spectrum is found to depend markedly on the concentration of sodium, 
Even at quite low concentrations (of the order of 0-001 gm/litre of hydrogen in a 
flame which consumes 1°5 litres/min. of hydrogen) the D lines are strongly reversed, 
and continuous emission extending for several hundreds of Angstréms on either 
side of the D lines is developed. At high concentrations, such as may be obtained 
by burning sodium metal in an atmosphere of air enriched with oxygen, there is 
strong emission from the continuum throughout the range 5500 to 7000a. 
(Crawford 1944, Hamada 1933). The visible radiation from all the yellow flames 
examined consisted primarily of emission from the D lines and from the associated 
continuum. 


Red flames. 


Strontium compounds are used to colour the red flames, but there are still some 
obscure features in the spectra. ‘The simultaneous presence of strontium and 
chlorine compounds in the compositions leads to strong emission of the red system 
of SrCl (the violet system also appears weakly). However, there are also present 
two closely-spaced sequences degraded to the violet which do not belong to the 
SrCl system. The long wavelength edges of these sequences are at 6884°5 and 
6114-2a. In strontium flames containing no chlorine compound these heads are 
more strongly developed, and, in addition, the region between them can be seen to” 
be filled with a very close pattern of what appear to be rotational structure lines. 
‘The most likely emitter would seem to be SrO, and it is some support for this 
suggestion that a heavy-current positive-column discharge through SrO contained 
in a silica tube leads to the production of the same bands. 


White flames. 


Some of the possible ways of making ‘“‘ white’ flames will be obvious from what 
has gone before. ‘Three methods which have been used are: (i) to develop an 
extensive continuum, e.g. the sodium continuum, in emission, (ii) to excite an 
extensive discrete band-system, and (iii) to excite simultaneously two more or less 
complementary band-systems. ‘The most suitable system for (ii) appears to be 
that of BaO, which may be produced conveniently by the oxidation of aluminium 
or magnesium by barium nitrate, while for (iii) suitable blending of the emission 


from SrCl (red), CaCl (yellow) and BaCl (green) has been found to yield an almost 
white light. 


§3. COLOUR AND SPECTRAL DISTRIBUTION OF ENERGY 


The spectral distribution of energy in the light from these flames has not been 
quantitatively determined, but the effects of light from secondary emitters on the 
colours of the flames can be considered in general terms. The red flames depend- 
ing on the emission from SrCl and SrO emit also yellow light (due to Na), green 
(due to BaCl) and blue (due to Sr). Green flames may emit, besides the green light 
from BaCl, some yellow from Na, red from SrCl and nearly white light from BaO. 
Blue flames, which rely on the D and E band systems of CuCl, emit also green light 
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(from the A system of CuCl), yellow (from Na) and red (from K). Calculations 
have been made on the effects of “diluting” light of each of the relevant colours 
with light of the wavel:ngths due to secondary emitters. Actually seven 
wavelengths have been considered and computations made for mixtures of light 
of each pair of wavelengths in proportions covering the whole scale from 1:0 to 


0:1. These seven wavelengths cover both the main and the secondary emitters 
(see Table 1). 


able=t 
Colour Wavelength a. Corresponding emitters with wavelength in a. 
Red I 6900 | SrCl, ca. 6200-6800; SrO, ca. 6000-6900; 
Red II 6200 f CaCl, 5800-6380; K, 7699 and 7665. 
Yellow 5890 Na, 5890, 5896. 
Green I 5300 BaCl, 5050-5350. 
Green II 5100 f ; CuCl, A system, ca. 5100-5500. 
Blue I 4600 | : 
Blue II 4100 f Sr, 4607; CuCl, D and E systems, 4200-4600. 


The method of computation is as follows. From the C.I.E. trichromatic 
coefficients for spectral light of given wavelengths, the corresponding specification 
for a series of mixtures is found (Wright 1944, Moon 1938). Let the trichromatic 
coefficients of two spectral colours to be mixed be (x,y, 21), (%2, V2, %2), and those 
of the resulting mixture be (x, Vm,%,).. Then if «, energy units of the first 
spectral colour be mixed with «, energy units of the second (where «, +a, =1), the 
trichromatic coefficients of the mixture will be given by x, =a 1%, +a %, 
Din = %1V1 + Xa Vo- 

Having found the trichromatic coefficients for a given mixture in this way, 
one finds the dominant wavelength with the aid of the C.L.E. diagram. A straight 
line is drawn from the point on the diagram representing equal-energy white light 
(x,, =0-333, y,, =0-333) to the point representing the mixture (Xin, Vm) and produced 
to meet the locus of spectral colours (or, for purples, the line joining the extremities 
ofthe locus). Let the point of intersection be («,,y,). The dominant wavelength 
of the mixture is given by the value of A corresponding to (x,, y,), the relative 
luminous intensity, by the value of y,,, and the colorimetric purity by the formulae 


(Hardy 1936) 
(Ps UY) Syria: Dw ge 0 aS ag 
Ym Vs —Vw Vm %s3— Ny 


The dominant wavelength, colorimetric purity, and relative intensity corre- 
spond roughly with the sensation attributes hue, saturation, and brilliancy, and 
give a fair idea of the appearance of a coloured light. The results are given in 
Table 2. 

Certain general deductions can be made from these figures, even though 
quantitative spectrophotometric data are wanting. (‘The variations of the several 
characteristics with the proportions of two colours are best seen by plotting them on 
graphs). As regards hue, the colours most easily spoiled by addition of light of 
other colours seem to be red and blue; as regards saturation, blue is much the most 
easily spoiled; as regards brilliancy, yellow and green. Blue light also causes the 
greatest decrease in purity when added to light of other colours. Addition of 
white light will reduce purity without altering dominant wavelength ; the emission 
of the BaO spectrum will have some such effect on green flames. 

3 ce 
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Table 2 


100 80 60 40 20 0 


100. 80 60 40 20 


Colour pair 


Red I (6900 a.) —yellow (5890 a.) 


Red II (6200 a.) —yellow (5890 a.) 


d (a) 6900 6256 6110 6018 5949 5890 | 6200 6104 6036 5980 5932 5890 
P 100 100 100 100 100 #8100 100 100 100° 100~ 1002 (am 
ym 2656 2985 3314 3643. 3972 4301 } 3083 3326 3570 3814 4058 43055 
Colour pair Red I (6900 a.)—green I (5300 a.) Red II (6200 a.)—green I (5300 a.) i 
d(a.) 6900 5996 5797 5627 5450 5300 | 6200 5926 5756 5598 5439 5300 
P 100 98:0 97:0 96:7 96:8 100 | 100 98-2 96:8 96:4 97:1 109 
ym 2656 3737 4317 5897 6978 8059 | 3083 4078 5074 6069 7064 8059) 
Colour pair Red I (6900 a.)—green II (5100 a.) Red IE (6200 a.)—green II (5100 4.) 
d(A.) 6900 6010 5767 5471 5187 5100 | 6200 5930 5715 5430 5177 5100: 
P 100 85-4 79-7 75:5 83-4 100 100 88-1 80-7 765 83-1 .10%R 
ym 2656 3625 4595 5564 6533 7502 | 3083 3967 4851 5734 6618 7502 
Colour pair Red I (6900 a.) —blue I (4600 a.) Red II (6200 a.) —blue I (4600 a.) 
d(a,) 6900 *4969 *5049 *5412 *5645 4600 | 6200 *4958 *5050 *5460 *5653 4600 
P 100 95:0 87-8 75-4 21-8 100 100 70:2 647° 55-0 31:8 106 
ym 2656 2184 1713° 1241 769 297 | 3083 2526 1969 1412 855 20m 
Colour pair Red I (6900 a.)—blue II (4100 a.) Red II (6200 a.) —blue II (4100 a.) ; 
d (a) 6900 *4970 *5050 *5370 *5613 4100 | 6200 *4960 *5048 *5420 *5626 4100 
P 100 100 100 100 98 +100 100 72:9 73:6 72-4 73:0 100 
ym 2656 2135 1613 1091 569 48 | 3083 2476 1869 1262 655 48 
Colour pair Yellow (5890 a.)—green I (5300 a.) Yellow (5890 a.)— green II (5100a.) 
(A) 5890 5766 5647 5525 5403 5300 | 5890 5749 5571 5323 5162 5100 
P 100 98:5 97:5 97:0 97:6 100 100 90-4 85-4 80:0 86-7 100 
ym 4301 5053 5804 6555 7307 8059 | 4301 5209 6444 7951 8205 7502 
Colour pair Yellow (5890 a.) — blue I (4600 a.) Yellow (5890 a.) —blue IT (4100 a.) 
d(A.) 5890 6013 *5048 *5598 *5693 4600 | 5890 6037 *5045 *5548 *5662 4100 
P 100,052 -7) 102222. 2017-20 Bez) S00 100: 52:2, ~ 27-0"  27-9° —28-Seae 
ym 4301 3500 2700 1899 1098 297 | 4301 3451 2600 1749. 898 


48, 


Colour pair 
A (A.) 
iP 


ym 


Green I (5300 a.) —blue I (4600 a.) 


5300 5172 5037 4922 4806 4600 
100. 84:4 70:0 63:1 57-8 100 
8059 6506 4953 3402 1850 297 


—_ 


Green I (5300 a.)—blue II (4100 a.) . 


5300 5178 5037 4911 4765 4100: 
100 S108 967-0) 956-25 42-9 100: 
8059 6457 4854 8252 1650 48. 


Colour pair 


Green IT (5100 a.) —blue I (4600 a.) 


Green II (5100 a.) — blue II (4100 a.) 


Ad (A.) 5100 5043 4980 4909 4814 4600 | 5100 5043 4977 4900 4781 4100 

P 100 94:0 88:0 81:8 74-6 100 100. 93-2 585-65 75°85 959 100 

ym 7502 6061 4620 3179 1738 297 | 7502 6394 4807 2950 1130 48 
Colour pair Green I (5300 a.) — green II (5100 a.) Blue I (4600 a.) — blue IT (4100 a.) 

A (A.) 5300 5247 5197 5158 5126 5100 | 4600 4560 4510 4450 4340 4100 

P 106 «97-3 ° 95-8 96:0 97°5 100 100 96 94 91 83 100 

vm 8059 8271 8333 8174 7867 7502 297 239 196 138 82 48 


A=dominant wavelength; P=colorimetric purity %; ym=relative luminous intensity. 


* These mixtures are purples, and the wavelen 
complementary colours, 


gths marked with asterisks are those of the 
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Red light is affected progressively less, both in hue and in purity, by addition of 
blue, green and yellow light. Yellow light j is altered in hue most by red or blue 
light; ; its purity is reduced considerably by addition of blue light but not appreci- 
ably by red. Green light is not greatly altered in hue by addition of up to 30% 
of light of any colour. Its purity is considerably reduced, more by blue than by 
red or yellow light. (On mixing green I and green II, the purity drops to about 
95 % asaminimum.) Blue light is altered in hue most by red, and in purity most 
by yellow light; the relative intensity is greatly raised by Badition of light of any 
of the other colours. The applicability of these general conclusions to the 
individual spectra is obvious. 

§4. CONCLUSIONS 

To be satisfactory, light from coloured flames should fulfil three conditions: 
(i) it should be of distinctive spectral hue; (ii) it should be saturated ; (iii) it should 
be intense. 

Thus the intense emission of the characteristic radiation of some molecular or 
atomic species which happen to be confined to a narrow range of desired wave- 
lengths would constitute a satisfactory solution. A relevant consideration here is 
the mode of excitation of the emitting species. On this we have as yet no direct 
information, but it seems most probable that the excitation in these flames is often 
purely thermal. If this is so, we can immediately deduce limiting relations for 
the variation of energy flux with wavelength and with temperature. For, quite 
generally, the energy radiated in unit time per atom in state 2 is A” hy,,,, where 
A” is Einstein’s coefficient of spontaneous emission for the change from state n to 
‘state m, which results in the appearance of a spectral line frequency »,,,. Now 
the values of A” for fully allowed transitions are not likely to vary much (cf. 

Mitchell and Zemansky 1934, Gaydon 1942); for a series, then, of ideal mono- 
chromatic emitters of different wavelengths maintained in thermal equilibrium, 
the intensities will be proportional to v,,,, exp (— hv/RT), where the exponential is 
simply the Boltzmann factor determining the population of the initial level 7 of the 
‘transition. To obtain relative luminous intensities we multiply this expression by 

the relative luminosity factor (Judd 1931). 


Table 3. Relative Luminous Intensities of Ideal Monochromatic Emitters 
for Thermal Excitation 


Colour Red Orange Yellow Green Blue Violet 
A (A.) 6800 6200 5800 5300 4500 4000 
POON nomemmtey << 106%) 5-252105" 2-6iXd0-" 2-710 1:1 x 105 2°35 x 105 
ae 00mm 45410 mom SO Om" cor? <dOs = 2-1 x 1052 69-7 105? | 1-o x 10s? 
i= 3 ())0> i. 0-22 eT 3°9 1:9 2: (<a Vie OO SE all Ogee 


A selection of values is givenin Table 3. The big variation with temperature 
of the intensity at any given wavelength was of course to be expected. More 
interesting are the facts that ideal red, orange, yellow and green sources have nearly 
the same relative luminous intensities at constant temperature, and that the blue 
sources are roughly 100 to 10,000 times less luminous. Any hopes of producing 
a blue flame depending on thermal excitation comparable in intensity with the 
other coloured flames are therefore illusory. 
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Quite generally, however, improvements are to be sought in cwo directions— t 
in increasing the intensity and in improving the quality of the light. The former 
implies an increase in the number of emitting species in their upper electronic 
state, which will generally be brought about by increasing the temperature of the 
flame. The effect of temperature may, however, be complicated if the emitters. _ 
are themselves formed by reactions taking place in the flame, e.g. BaCl,>BaCl+ Cl; 
or if decomposition becomes important at high temperature, e.g. BaCl> Ba+ Cl; 
or if competing equilibria are important, e.g. BaCl+O—-BaO+Cl. For a given 
emitter, the quality of the colour can be improved by reducing the light-output 
from undesirable secondary emitters. 

Of existing compositions designed to give coloured flames, those giving red 
flames (based on SrCl and SrO) are the most satisfactory, both from the point of — 
view of intensity and of colour. This is readily understandable in terms of the © 
above discussion. It so happens that only the red systems of SrO and SrCl are of 
much importance; the blue and near ultraviolet systems of both molecules are — 
developed (as might be expected) only very weakly. Furthermore, the only — 
important secondary emission is that from the resonance lines of potassium which ~ 
comes from the oxidizing agents generally used—potassium chlorate or perchlorate. 
These lines lie in the red region of the spectrum and, therefore, contribute to the — 
desired colour. The sodium D lines appear but do not seriously affect the hue or © 
saturation (cf. Table 2). The best results with this type of composition cannot, 
however, be obtained without purification of the components from barium. 

A sodium source is the only likely solution to the problem of finding a yellow — 
flame, and some observations suggest that a similar source may also give a nearly ~ 
white light of high intensity. It is not difficult to reduce the effects of red and blue ~ 
light from secondary emitters to negligible proportions. 

The distinctive colours which remain are green and blue. Intense flames of 
these colours do not appear to have been produced. It has been shown above that 
if the excitation is thermal, the operation of the Boltzmann factor and of the 
luminosity function combine to reduce the theoretical intensity of a blue flame to. 
quite a small fraction of that expected for, say, an orange flame. Moreover the 
saturation of a blue flame is much reduced by relatively small proportions of light 
from secondary emitters ; the quality of the colour from existing blue compositions. 
is adversely affected by weaker emission from the A, B and C systems of CuCl, and 
perhaps by strong emission of the K lines, which, however, have a wavelength near 
the limit of visibility. Unfortunately it is not easy to find satisfactory alternatives. 
either to the potassium salts or to CuCl as emitter. 

The most serious difficulty encountered in attempts to improve the green flames 
based on BaCl is the emission from the extensive system of BaO. (Small amounts. - 
of Na, Ca and Sr also affect the colour adversely, but these could largely be 
removed by careful purification.) Our results suggest that the ratio BaO/BaCl 
is greater the more intense the source. Simultaneous improvements in quality 
of colour and in intensity are not therefore to be expected in sources which rely on 
emission from BaCl in an atmosphere which contains free oxygen or oxidizing 
agent. Somewhat surprisingly, there seem to be no very promising alternatives 
to BaCl among the diatomic molecules, although BO might be useful: it is also 


possible that a suitable means might be found for exciting the green line of thallium 
(5350 .). 


# 


Ee a 
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_ ABSTRACT. From an investigation of the discharge spread in Geiger counters filled 
_ with hydrogen, neon, argon or helium, it is possible to show that in certain circumstances 
the discharge spread is due almost entirely to photo-ionization in the gaseous counter 
filling. The results enable the absorption coefficients of the operative radiations to be 
found. 


Sie. DN RODUETION 


certain forms of electrical discharges in gases (Meek 1940, Loeb and Meek 

1941, Weissler 1943). It is also of great importance in upper atmosphere 
problems (Sayers 1943, Massey 1938, and many papers by Appleton and his 
collaborators) and finally is relevant to work on Geiger counter discharges (Korff 
1946, and many others). 

Price (1943) has pointed out the singular lack of experimental data on photo- 
ionization in the elementary gases, a situation which is in contrast with that 
obtaining in studies of alkali metal vapours (Ditchburn 1928 etc., Lawrence and 
collaborators, and others). Page (1939) has published an admirable review (with 
references) of the theoretical and experimental work up to that time, with particular 


Pp HOTO-IONIZATION in gases is of interest, particularly in relation to studies of 
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: 
reference to work on alkalis (Cs, K, etc.), but omitting reference to certain important | 
theoretical papers (Wheeler 1933, Vinti 1933). The later theoretical work is due 
largely to Massey and his collaborators, particularly Bates (1939, 1947 etc.). 

It was thought therefore that further experiments on the elementary gases, 
employing a variety of widely different techniques, might be of value and because 
of certain advantages of the method, it was decided first to investigate and refine 
a Geiger counter method, due originally in a simple form to Greiner (1933) and 
recently used in experiments on complex gas mixtures (CH,, C,H;OH/argon etc.) 
by Craggs and Jaffe (1947 a, b), Balakrishnan, Craggs and Jaffe (1948), Alder and 
his collaborators (1947) and Liebson (1947). The above workers used different iL: 
techniques as may be seen from reference to the papers. 

The method is shown to suffer from certain important limitations which are | 
discussed in §6, but the present Geiger counter experiments differ from previous __ 
comparable work in that photo-emission from the tube walls is eliminated as 1 
described below in §3. The wall effect which was present in Greiner’s (1933) 
and Christoph’s (1937) experiments might, it was thought, have exerted an ~ 
influence and rendered the results spurious. 


§2. SUMMARY OF PREVIOUS THEORETICAL AND 
EXPERIMENTAL WORK i 

Page’s review (1939) discusses in some detail the interest evoked in the subject 
by Kramers’ (1923) theory, originally worked out for hydrogenic atoms, and the 
limited extent to which it agrees with experiments for the alkali metals. Trumpy’s 
later (1931) wave mechanical treatment of the latter was more successful (Page — 
1939) and gives reasonable agreement with experimental data provided, for — 
example, by Ditchburn and his school (Ditchburn 1928, Braddick and Ditchburn . 
1934). 

Wheeler (1933) and Vinti (1933) provided wave mechanical treatments for 
helium. Vinti shows, for example, that the atomic cross-section for photo- 
ionization at the principal series limit (A=504 a.) is about 8-4 x 10-18 cm?, falling . 
off to about one-tenth of this value at about 180a. The cross-section for the 
carbon Kz« radiation (44-6 a.) was measured by Dershem and Schein (1931) and 
agrees satisfactorily with Vinti’s calculated value. 

The atomic cross-section o, is related simply to the linear absorption coefficient 
(pressure p) by o, =/No, where No is the number of atoms per cm? (pressure fp), 
and pis given by J=I,e~”, where J, and J are the incident and emergent intensities 
for a thickness x of absorber. 

Theoretical treatments for atomic hydrogen are due to Suguira (1927), 
Gaunt (1929), Menzel -and Pekeris (1935) and others. Massey (1938) shows 
o, 0-6 x 10-1’ cm? at the series Jimit. 

Since the writing of Page’s review, Bates’ paper (1939) appeared in which 
results of great interest are described. ‘The continuous absorption coefficients 
per 2p electron for B, C, N, O, F, and Ne are calculated for various ejected electron 
energies. ‘lhe introduction of the appropriate weighting factors gives the 


following values of absorption cross-sections at spectral heads. (Bates, private 
communication). 


w 4° 


Table 1 


Atom B Cc N O F Ne 
Cross-section (107!? cm?) 1:9 2:4 ED) 0:45 0:61 0:58 


Some Experiments on Photo-ionization in Gases 4I 


It will be seen that Ne is the only convenient element with which the above 
| results may be checked. Bates’ results show that the Ne absorption varies much 
_ less with energy than does the helium absorption worked out by Wheeler (1933) and 

Vinti (1933). Experiments to check these He and Ne results are in progress at the 
| University of Liverpool, using a different method from that to be described below. 
| Experimental data on the elementary gases are scarce. Greiner (1933) and 
later Christoph (1937) used a Geiger counter method in which the discharge of 
one counter (a) was allowed to initiate a discharge in another counter (b) enclosed 
' in the same chamber, by virtue of photo-ionization produced in (b) by energetic 
_ photons from the discharge in (a). Their results are summarized in Table 6. 
_ Raether (1938) described a novel method in which he measured photo-ionization 
caused in the gaseous filling of a Wilson cloud chamber, by photons passing 
through a very thin cellophane window from an external corona discharge into the 
_ chamber. Geballe (1944) measured the absorption coefficient for highly absorb- 
able radiations produced in a Townsend discharge and detected photo-electrically 
at brass surfaces placed at different distances from the discharge. These results 
are also given in Table 6 and the discussion in §6. Finally, Schneider (1940) 
_ measured absorption coefficients in air for wavelengths lying between 380 and 
1600 ., using a vacuum spectrometer. 


§3. APPARATUS 
A new type of divided cathode, single anode counter—with provision for 
collimating the photons originating in the discharge—was constructed (Craggs 
and Jaffe 1947a, b). It consists of six copper cathodes, each 5cm. long, 2cm. 
internal diameter and spaced 1-2cm. apart, with an anode wire of 5 mil tungsten 
_ (Figure 1). 

In order to eliminate the cathode effect by collimating the photon beam to the 
neighbourhood of the wire, a window of C40 glass (which seals well to Fe/Ni/Co 
alloy), was attached to the front of each except the first cathode. The central hole 
in each window was 3 mm. in diameter in the beginning. Later, that for cylinder 
B was replaced with another window with a hole of diameter 1mm. only. The 
alignment of the five windows was carried out with great care and the collimation 
checked before the wire was inserted and the tube sealed at the ends. 

A special counter was used in order to check that no radiation passed through 
the glass of the windows used in the previous counter. This consists of two separ- 
ate compartments separated by a circular piece of C40 glass of the same thickness 
as that used in the six-cylinder window counter. A cylindrical piece of Kovar wire 
is fused to the glass through a hole in the centre. ‘The copper cathodes (above 
dimensions) are spaced 10cm. apart. 

The quenching circuit used, shown in Figure 2, is similar to that of Neher and 
Harper (Craggs and Jaffe 1947b). Careful tests were made to ensure that 
the coincidence circuits were symmetrical and that both channels were identical in 
sensitivity. The two-channel quenching circuit, which was arranged to take 
pulses from the cathodes of the divided counter, supplied the pulses to a scaling 
circuit and mechanical counter. Gases of cylinder grade were used.* 


* The experiments were intended to give information relevant to the conditions of sparks in 
cylinder gases, but it is necessary to stress that the observed values (Table 6) could be given 
by absorption in impurities (having p’~40cm™! at 10 cm. Hg) being present to the extent of 
~0:25%. The Hz and Ne results should be affected less than those with A since in the former 
cases the main impurities have higher ionization potentials than the main gas (see also Raether 


1938). 
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§4. PRELIMINARY TESTS 


It was necessary to determine the resolving time of the coincidence circuit used, 
in order to check that it was short enough to reduce random coincidences to a 
reasonably low level, but was not comparable with or less than, the time for dis- 
charge spread, which would be of the order of 10 sec. for photon induced spread. 
Two similar separate counters, each equivalent to one section of the counter of — 
Figure 1, were set up 10cm. apart. By taking the individual counts N, and N,, 
when one of them is irradiated with a collimated source of gamma rays from radium 
and the chance coincidence counts C in the same conditions, the resolving time, 7, 
can be calculated from the equation C=2N,N,7 + K, where K, the true cosmic ray 
coincidence count, may be eliminated by using two sets of data. N, and N, 
included the counter background in each case. By reducing the operating voltage _ 
of one counter just below its starting potential, the absence of circuit pick-up 
effects was checked. 

The following results were obtained in a typical experiment with two separate 
counters : 


(a) With no gamma radiation: N, =648/min., N,=549/min., C =2-8/min. 
(b) With gamma radiation: N,=906/min., N,=1321/min., C=6-9/min. 


7 was calculated and found to be about 1:5 x 10~*sec. 

The time constant involved in the circuit of Figure 2 is such that coincidences. 
are missed at high counting rates, and the effect was shown by calibration curves 
relating to the number of counts per cylinder and the coincidence counts. ‘To 
avoid distortion of the discharge spread curves (Figures 4-6) by such a spurious 
effect, the runs for the latter figures were all made at as constant a counting rate as 
possible. 

Experiments were conducted with the test counter to check that the glass used 
for the windows is the six-cylinder counter was opaque to the radiation originating 
in the gas discharges. ‘The two-chamber counter was filled with argon and hydro- 
gen in turn at various pressures and one of the cylinders was irradiated with a 
collimated source of gamma rays. Determinations were made of the counting 
rates of this irradiated cylinder and of the other in the neighbouring compartment 
together with the coincidence counts for the two cylinders. These rates, when 
compared with the corresponding rates in the case of two individual counters 
arranged at the same distance apart, show that the coincidence rate for the two- 
chamber counter is equal to the chance coincidence for the two separate counters. 
Hence it may be safely assumed that the glass used for the windows in the six- 
cylinder counter (F_gure 1) is opaque to the radiation originating in the discharge 
and so the photon beam is not irradiating the cathodes. Table 2 gives a summary 
of the data observed. ‘The counting rates for two cylinders (in the six-cylinder 
counter) whose distance of separation is approximately equal to that of the cylinders 
in the two-chamber counters are also recorded in this table. This shows the 
amount of discharge spread at the distance, i.e. at approximately 10 cm. 

Previous work (Craggs and Jaffe 1947b and references cited there to earlier 
experiments by other workers) showed that the discharge spread in H, and A for 
divided cathode counters of the dimensions shown in Figure 1, but without colli- 
mating windows, was approximately 100%. Greiner (1933) who did not employ 
collimation found less than 100°% spreading, but this discrepancy is explicable by 
different cathode work functions in Greiner’s and the present experiments. 
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§5. EXPERIMENTS WITH SIX-CYLINDER WINDOW COUNTER 


The experimental arrangement is shown diagrammatically in Figure 3. The 
counter was filled, at a definite pressure, with the gas under investigation and’ 


Table 2 
Type of Gas Pressure 
counter filling (cm. Hg) (a) (b) (¢) 
‘Two separate Argon 10 121 906 6-9 
counters Argon 5 1681 1496 21 
‘Two-chamber Argon 10 638 594 By) 
counter Hydrogen 10 640 ' 415 8 
Hydrogen 5 1209 80 3 
Six-cylinder Argon 20 419 208 87 
window counter . 10 1015 1071 661 
(Figure 1) 5 807 754 607 
2-5 600 494 384 
Hydrogen 20 443 152 135 
| 10 466 323° 307 
| 5 382 350 332 
25 361 360 331 


(a) Counting rate/min. for cylinder irradiated with gamma rays. 
(b) Counting rate/min. for other cylinder. 
(c) Coincidence counting rate/min. 


O-1cm.diahole O-3cm.diathole ; mall, 
; | hi stabilized supply <—__ _ 
NSS NK sie 
1 Nis \= ses is a 
a L = = A= =, = 
as a/\ ~O? © zm © \ © Output - 
O-lom thick glass \ Fe/Ni/Co 0-013cm dia 
window Alloy Cap tungsten wire 


200000 Mn 


Figure 1. 6 cylinder window counter. 


(eee 25cm, —>| 0-5 mgm | 


radium 


ai 


1 
| Window Counter 


200000 Ma 


To Quenching Circuits 


Figure 3. Diagram of experimental arrangement. Figure 2. T'wo-channel external 
quenching circuit. 


mounted before a lead collimator so that cylinder A was irradiated with a collimated’ 
beam of gamma rays from 0:5 mgm. of radium embedded in the lead collimator. 
Cylinders A and B were connected to the counting channels and their counting: 
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rates were taken. The coincidence rate for the two cylinders was also recorded. 
“The measurements were repeated with cylinder A and each of the successive cylin- 
-ders, C, D, E and F connected to the set. 

If N, and N, are the counting rates of two cylinders and Ng the coincidence 
‘counting rate, the total of the true counts of the two cylinders will be N, +N; — No, 
-since the spreading of a pulse from the first to the second will cause an apparent 
increase in the counting rate of the second and so also for the pulses which spread 
-from the second to the first. The percentage pulses which cause coincidence is 

100No/(N,+N2—No). 

The percentage spreading was calculated for each cylinder and plotted against 
‘the distance of that cylinder from the irradiated cylinder. 

In order to check the absence of electron diffusion and to ascertain if any 
spurious charges which might affect the results were accumulated on the inter- 
mediate cylinders when the remote ones were counting, counting rates were taken 
with the intermediate cylinders (a) floating, (b) at a potential less than the starting 
potential and (c) shorted to the wire. The data are recorded in Table 3 and it is 
-seen that the results are not affected in any way by leaving the intermediate cylinders 
floating.” 

The curves obtained were flat topped, showing that there was 100°, spreading 
for the first two or three cylinders. This may be due to an excess of photons 

-generated in cylinder A, since one photo-ionizing photon will give a count. 


‘Table 3 
Intermediate cylinders Floating At 700 v. Floating At500v. Floating pas 8 
Counting rate/min. of A 948 947 =P i35 i 374 170 183 
Counting rate/min. of B 942 958 358 352) 189 167 
Coincidence counting rate/min. 587 624 151 147 116 123 


The counter was now slightly modified. The window in cylinder B had a hole 
‘of diameter approximately 3mm. This was replaced with another, with a hole 
of diameter 1 mm., so that less photons might spread into the other cylinders, and 
measurements were repeated. Although the discharge spread curves (_.e. 
coincidence rate—cylinder number (Figure 1) with A irradiated) for argon were 
straight lines, those for hydrogen were flat topped with nearly 100°, spreading 
in the initial stages, showing that more than one ion-pair was produced for each 
pulse in cylinder A and so excluding any measurements of absorption coefficients. 

Another objection to these early experiments was the possible production of 
photo-emission from the edges of the glass windows, which could then trigger off 
discharges in the corresponding cylinders. If the windows were emitting 
photo-electrons, then an asymmetry in discharge spreading would be noticed, 
for example, cylinder A would give more discharges in cylinder C than C would 
givein A. 
The asymmetry in spreading was detected experimentally. Cylinder A was 
irradiated with the collimated source of gamma rays, and the counting rate, when 
this cylinder alone was connected to the set, was taken. Next cylinder E was also 
‘connected to the set, and the counting rate for A with E and E with A and the 

* Later experiments on the velocity of discharge propagation (to be published) indicate that 


-other diffusion phenomena are not important. One of us (J. D.C.) wishes to acknowledge a 
valuable private communication from Professor L. B. Loeb on these matters. 
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| coincidence rate for them were found. The collimator was then moved so that 
| the gamma ray beam irradiated cylinder E and the counting rate of E (with A 
| disconnected) was adjusted to be the same as that of A (with E disconnected), by 
| manipulating the position of the radium in the collimator. Then the counting rate 
| of E alone, of E with A connected to the set and the coincidence counting rate were 
| measured. Since the individual counting rates of A and E were adjusted to be 
' equal, if the spreading was symmetrical, the counting rate for A with E and E with 
A as well as the coincidence rates in the two cases of irradiation should be equal. 

‘Table 4 shows the data and it is easily seen that the spreading is asymmetrical. 
_ The degrees of collimation for the two directions of discharge spread should be. 
- sensibly equal. ; 


Table 4. Counting Rates for various Conditions 


A, with E E, with A A, with E E, with A Coincidence 


disconnected disconnected connected connected counts 
A 426 41 423 330 247 
E 33 453 87 513 69 


First column indicates cylinder irradiated. 
Figures give number of counts/minute. 
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Figure 4. Discharge spread Figure 5. Dischargespread Figure 7. Absorption coefficients 
in hydrogen. in argon. (arbitrary units) plotted as a 
function of pressure. 


In view of this asymmetry in spreading all subsequent experiments were 
modified to eliminate this window-edge effect. Instead of irradiating cylinder A 
_ and measuring the spreading in the other cylinders, each of the cylinders B, C,... F 
were irradiated successively with the collimated beam of gamma rays, and the 
spreading to cylinder A was measured in each case. ‘The spreading curves for 
hydrogen, argon and neon are plotted in Figures 4, 5 and 6 respectively, at various 
pressures; Figure 7 shows a plot of the slope of the spreading curves against pressure. 
The apparent absorption at zero pressure is presumably due to the imperfect 
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collimation. Table 5 gives a summary of the starting voltage and the working © 
voltage in each case. Experiments made in the initial stages show that the absorp- — 
tion does not vary with overvoltage although there is greater spreading at higher 
voltages. z 


‘Fable 5 
Gas p vg UW Gas p Us Uw ‘ 
Argon 98% 95} 580 620 Hydrogen 99% 2:5 720 760 
(~2% Ng) 5 760 800 (~1% 0.) 5 880 920 be, 
10 920 960 10 1160 1200 . 
20 1220 1260 20 1600 1640 ¢ 
-Neon 98 % 5 380 420 Helium 98% 10 1100 1200 | | 
(~2% He) 10 400 440 (~2% Ne) 3 830 930, : 
20 440 480 }) 


p=pressure (cm. Hg); vg=starting voltage; vy—=working voltage. 

The values of absorption coefficient jz (see §2 for definition) may be derived 
‘from Figures 4, 5 and 6 and are given (for 10cm. Hg pressure) in Table6; they were 
quite close to those obtained in the presence of emission from the edges of the glass 
collimating windows. Data for helium were also obtained in some preliminary 


experiments, but owing to lack of helium, further data without the window effect 


ve heer aie” ie, ae nea 


Table 6. Values of Absorption Coefficient » for p=10cm. Hg 
Experi- : 
H, Air O, H N O He Ne A p Author eer. 
remarks i 
0 ee ee eee — ~0-1 Schneider Vacuum : 
(1940) spectro- 
meter 
0-11 0:26 0-66 14-50 Raether Wilson 
(1939) chamber 
55 —- —- —- — ~~ — — _ 0-1 Geballe Photo-cell 
(1943) (brass 
cathode) 
018 — — ~ ~~ — — a ~10 Greiner Geiger 
(1933) counter 
Soy 70) —- —- ~— — _- — -1'55 Christoph Geiger 
(1937) counter 
—_—- — 22 -—- — — Massey Calculated 
(1938) at series 
limit 
_- — 80 16 — 21 a — Bates Calculated 
(1939) at series 
limit 
— ee — a — Vinti Calculated 
(1933) at series 
limit 
0:12 — ee ee () 5) 5) le) 25 — | maT eserat Geiger 
authors counter 


p=gas pressure (cm. Hg) used in experiments. 
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(see above) were not obtained. However, it seems reasonable to state that the 
_ absorption coefficients for helium lie between those of hydrogen and neon. 
In some experiments with a long counter (Figure 8) cylinder A was irradiated 


) with the collimated radium source and the discharge spread in cylinders D, E,...H 
7) was measured in the usual way. Owing to the long distances involved in this 


1 


® 0-3cm.dia. hole © @ eo) o) 


Fe/Ni/Co 
‘ ae Alloy Cap 
0-Icm. thick glass window 0-013cm. dia. tungsten wire 


Figure 8. Long collimating counter. 


counter, there was perceptible spread only in neon and hydrogen at low pressures. 
_ At higher pressures in these gases as well as in argon there was almost complete 
absorption and practically no spreading at such great distances. 


§6. DISCUSSION OF RESULTS 

Previous experimental workers have often tended to ignore the theoretical 
work on He (Vinti 1933) and atomic H (Table 6 and elsewhere) and arranged - 
experiments which obviously would not have permitted such radiations to be 
measured, due to the use of excessive path lengths, etc. For example, taking 
Pair 28 50 cm7! (Schneider 1940, Table 6), J)/7~100 for «=0-1 mm. at 1 atmo- 
sphere pressure and the absorption at 1 cm. Hgand x~1cm. would also be about 
100: 1, so that counter experiments could not be expected to yield much useful 


} information on such absorbable radiations. "The same conclusion could have been 


reached from the 1933 work on He and from Bates’ (1939) paper. These high 
values of absorption coefficient lend qualitative but strong support to the streamer 
theory of spark breakdown proposed by Loeb and Meek (1941 and references there 
cited). 

However, since the values of o, with frequency for He and Ne are different, it 
was hoped that Geiger counter experiments might demonstrate such differences, 
and indeed for this reason it was also thought that the neon values foro, might be 
useful because, since a gaseous discharge (Geiger counter or Townsend [Geballe] 
type) would generate photons with a finite width of spectrum, the effect of such a 
spectrum in determining the ‘‘mean”’ absorption coefficient of the radiations 
should be less in Ne than in other gases. 

If a suitable monochromatic source of photo-ionizing radiation could be found, 
and the x-rays used by Dershem and Schein were of too short a wavelength to give 
much valuable information on absorption near the series limit, then the filtering 
effect of distance would clearly be of less importance. 

It is concluded that the vacuum spectrometer method has exceptional value 
(Table 6) since low gas pressure (<1 mm. Hg) may be used for absorption, and most 
electrical discharge devices (e.g. Geiger counter) need a higher gas pressure than 
this, although Geballe’s Townsend discharge could be operated at approximately 
Imm. Hg. Geballe’s work might be criticized that the emissions from two 
(different) photo-electric surfaces of brass (not degassed) at different distances 
from the photon generating discharge were compared without direct comparison 
of the work functions of the two surfaces. 
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§7. CONCLUSIONS 


In conclusion, therefore, it appears that the present (Geiger counter) method 
and those of Raether, Greiner etc. cannot be expected to give true values of the 
absorption coefficient for photo-ionizing radiations, due to excessive filtration 
effects. The existence, in such discharge devices as Geiger counters, of radiations 
having p~0-1 cm at 10 cm. Hg pressure shows (from the work of Wheeler (1933) 
and Vinti (1933) for helium) that radiations having A~50 a. may be present in such ~ 
discharges.* It is unfortunate that photo-ionization cross-sections for the common 
gases H,, N, and O, have not, owing to certain theoretical difficulties, been worked — 
out. Data for H, N and O (atomic) are given in Table 6. 
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/ ABSTRACT. The paper describes a study of the excitation temperature in certain spark 
«| discharges of accurately known current characteristics. 'The excitation temperatures are 
| found by measuring intensity ratios for certain spectral lines where the relevant transition 
| probabilities are known, and are comparable with earlier work by Ornstein and his colla- 
‘6 borators on arc discharges. 


The excitation temperatures are, as would be expected, higher than those generally 


¢ found for arc conditions. 


Certain peculiarities relating to the evaporation of electrode metal are described. In 


+ particular the presence of discontinuous evaporation was noticed, and it was found that 


different electrode metals showed considerable variations in their behaviour in this respect. 


§1. INTRODUCTION 


HE study of phenomena occurring in spark gaps is of interest for several. 
reasons; there are many technical applications of spark discharges for 
various and widely differing purposes, such as the ignition of explosive 


| gaseous mixtures in internal combustion engines and the spectroscopic analysis. 


of many substances. The behaviour and function of the metallic electrodes 
used in these circumstances are clearly of considerable importance, and many 
aspects of the discharge-gas—electrode-metal interaction are little understood. 

The present experiments were undertaken in order to obtain preliminary 


? information on excitation conditions in spark discharges in air whose electrical 


characteristics were defined more clearly than in previous comparable work 


; (see §§2 and 6); in particular, to give values of the excitation temperature in the 


discharges, using spectral lines of Ba and Mg for which the transition probabilities. 


t are known from separate experiments with arc discharges. ‘The latter work is 


due largely to Ornstein and his collaborators. 
Extensions of this work in various directions, particularly to the case of more 
arc-like discharges, are taking place. 


§2,. SUMMARY OF PREVIOUS EXPERIMENTS 

Unlike the excitation conditions in pure arcs (e.g. Ornstein and Brinkman 1934) 
and even in glow discharges, the modes of electron-ion or electron—atom inter- 
actions in pure spark discharges have: been little studied and are not so well 
understood, probably because of the relatively greater complexity of the discharge. 
It is difficult, in the first place, to define a pure spark but for present purposes the 
expression is taken to mean a transitory gaseous discharge in which (a) the voltage 
gradient is high and is changing with time, (4) thermal equilibrium is approached 
to an unknown, but probably limited, extent, (c) the general electrical state can be 
described as being between that of some forms of corona (streamers etc. (Loeb 
and Meek 1940)) and glows or arcs. The general physical characteristics of 
sparks are well known, viz. brightness, tortuous nature etc. 
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Previous work on such discharges, with particular relevance to the study of © 
spectroscopic analysis, is due to Kaiser and Wallraff (1939), Langstroth and 
McRae (1938 a, b) and many others. A useful and accurate review of the subject _ 
is by Llewellyn Jones (1945). In the above researches, oscillatory discharges have 
been used and photographic records of spectra generally give an integrated record 
of the emission for the whole wave train. This may be misleading, since the 
spark discharge probably only obtains for the first half-oscillation and tends to 
become an arc-like discharge during the latter oscillations of the train. For 
example, Langstroth and McRae (1938 a) used a circuit which gave a damped train 
of waves with a period of 100sec., and found from observations of the well- 
known cyanogen bands, and from certain tin lines (2850, 3009, 3262, 3175.) 
that thermal equilibrium was established during the discharge time, 1.e. the 
excitation was thermal with approximately equal gas, electron and ion temper- 
atures. The temperature found was about 9,500°K. However, it was pointed out 
(Langstroth and McRae 1938 b) that the potential difference across the spark gap, 
for every half-oscillation after the first, was about 300 volts falling rapidly 1 to about 
35 volts in a fraction of a period. Hence it is not surprising, in the present 
authors’ opinion, that thermal equilibrium was established. 

A few measurements (Langstroth and McRae 1938a) were made with the 
same circuit modified to give a damped wave train with a period of about 1 psec. | 
It was found that thermal equilibrium was not established since temperatures 
‘deduced from the tin lines lay between 8,700° and 11,000°x. Apparently no” 
measurements with the cyanogen bands were made. However, from these 
measurements the departure from thermal equilibrium (see § 6 for a discussion of 
energy distributions and statistical equilibria in certain high current density — 
-discharges) is apparently not very large, considering the difficulty of interpreting 
the data and the uncertain knowledge of the exact form of the discharge. This 
may again be due to the nature of the discharge, and the finite length of the oscilla- 
tory wave train. 

Recurrent sparks were used carrying square waves of current (Craggs, Haine 
and Meek 1947, Craggs and Meek 1946). In this case the wave form of the spark 
current is accurately known. About 20-30kv. (accurately measurable) is applied 
to a spark gap arranged in series with a special 3-electrode trigger gap (Craggs, 
Haine and Meek 1947) so that the test gap breaks down in much less than 1 psec. 
"The following current wave is sensibly rectangular (see Figure 4(e)) and of con- 
trollable duration, e.g. 1-9 usec. in the present experiments (§3). The recurrence 
rate is controllable and is usually about 50/sec. The voltage across the gap falls 
from its initial value (say 20 kv.) to about 250 volts (in air) in about 0-5 usec. 
(Craggs, Haine and Meek 1947, Figure 27). The sparks were observed with 
-electron multipher tubes, coupled to a cathode-ray oscillograph, and with aq 
rotating mirror camera. 

Further work will be carried out with even shorter discharge times and a 
separate detailed investigation of the voltage characteristics of short sparks is in 
progress. 


§3. DESCRIPTION OF THE APPARATUS 


The square current pulses were produced by the charging and discharging of 
an artificial transmission line using a circuit described by Craggs, Haine and Meek 
(1947). ‘The oscillograph and spark firing circuits were electrically synchronized. 


Short Sparks in Air ~ 51 


The optical system for spark observation was arranged so that light from the 
‘spark was focused on the collimator lens. ‘The telescope eyepiece was replaced by 
another slit in the focal plane of the spectral lines; the electron multiplier, which 
‘vas coupled through an amplifier to the cathode-ray oscillograph, was mounted 
just beyond the telescope slit, and the telescope slit width was adjusted to give 


dequate resolution. Further details of the amplifier are given by Craggs and 


The rotating mirror used for observation of the vapour clouds emitted by the 
“spark electrodes was such that a resolution in time of 40 ywsec/cm. was obtained on 
ithe photographic plate. ‘The plates were developed and microphotometered with 
; care using normal techniques, such as continuously brushing with cotton wool 
iduring development. 


§4. EXPERIMENTAL RESULTS 

In order to determine excitation temperatures by the measurement of spectral 
‘lines the transition probabilities of the latter must be known. Ornstein and 
+ Brinkman (1934) and others have measured the transition probabilities for many 
‘lines produced in arcs where, as can be shown by observation of the CN bands, 
thermal equilibrium is obtained and the temperature may be found. The 
} excitation temperature is found from the CN bands or in other ways, and this 
‘knowledge enables the relative transition probabilities of lines, whose relative 
intensities are known, to be found. The operative relation is 
T/T, =(AyGy1/A,G@ov2) exp ((E,—Ey)/RT}, «ss (1) 
| J, and J, are the measured intensities of two lines of frequencies v, and v, and 
/ relative transition probabilities A, and A,, G, and G, are the 1espective statistical 
) weights, E, and E, are the energy values (ergs) of the initial levels for the two 
‘transitions, and k and T are Boltzmann’s constant and the gas temperature 
respectively. We shall show later that even when thermal equilibrium is not 
| attained, T may be almost equal to 7, (electron temperature) if certain conditions 
| are satisfied by the discharge. 
The main purpose of the investigation was to derive information on spark 
-channels in air between metallic electrodes, using the known transition prob- 
| abilities, found from the work on the arc discharge. The present study is limited 
| to visible radiation; Ba and Mg seemed to be the most suitable metals of those 
whose transition probabilities have been measured. ‘The work on Ba and Mg 
, was carried out by Mason (1938), Kruithof (1943) and Kersten and Ornstein (1941). 
Tests on other metals, e.g. Ca and Cd, whose transition probabilities were known, 
showed that they were not so suitable for the present preliminary study. 
The electron multiplier experiments are first described. ‘The spectral 
' response (Figure 1) of the multiplier was such that many red lines (e.g. some of those 
of Bal) could not be measured. The Mg triplets 3838, 3832 and 3829a. and 
5184, 5173 and 5167a. were measured as groups of three added lines since the 

available dispersion was too small to give adequate resolution into the separate 
component lines. Figure 2 shows some typical results with these Mg triplets 

(approximately 100amp. pulse for 4:5 sec. with a 3-mm. gap in air between 
pointed electrodes). 

When tests with Ba metal electrodes were made, the spark-to-spark fluctuations 
in light output, for constant current pulse shape, were too great to allow satis- 
factory oscillograms to be taken. 'To find the source of these fluctuations, the 
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spark gap was moved across the optical axis of the system and the peak heights of 
the light-emission-time oscillograms (e.g. Figure 4(d)) noted for each of a number 
of fixed spark conditions. From the plot of peak light intensity against displace- 
ment from the axis it was possible to show that the effect of random spark movement 
(with fixed electrodes on the optical axis) on the oscillogram variations was small, 
The observed spark light fluctuations corresponded to off-axis movements of 
about 8mm. in a typical case (Figure 4(d)) whilst, in practice, the maximum 
random spark movements were only about 2mm. from the mean (axial) position. 
Much greater fluctuations than those of Figure 4(d) were often observed. It 
seems therefore that the great observed fluctuations in the amount of light (at 
any particular wavelength) emitted in a single spark of the present character must 
be due to electrode effects. Ba serves as an extreme case, where the fluctuations 
are most pronounced. It must be emphasized again that the current wave from 
spark to spark showed no fluctuations. 


x 10 April 1947 
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Figure 2. Light/time and current/time — 
oscillograms taken with Mg electrodes — 
and 4 usec. current pulse. 
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temperatures are derived from the ratio. 
of intensities of the MglI triplets. 
A= 3838 a. and A=5184 a. 


It is necessary to show that the light emission fluctuations as in Figure 4 (d) are 
not due to multiplier noise. The relative noise fluctuations (R,) are given by 
(Dieke et al. 1946) R,=1/,/(nt), where ¢ is the time interval considered, and n is 
the number of electrons flowing per second. The multiplier photo-cathode 
current is not less than 20x 10-” amp. and for t=10-®sec., R,~9%. The 
fluctuations shown in Figure 4(d) (~30%) are much greater than this and are in 
fact themselves much less than the fluctuations observed in other cases. Further 
confirmation of the relative unimportance of multiplier noise is provided by 
observations on hydrogen sparks which show fluctuations considerably less than 
9% in the conditions of the present experiments. 

‘The use of various barium salts in drilled electrodes gave no improvement, but 


specially prepared alloys Ba/Cd (20%/80%) and Ba/Mg (20°%/80%), gave 
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“successful results. The light-time oscillograms still showed fairly large fluctua- 
sions, but enabled satisfactory readings to be taken. Fitzgerald and Sawyer 
1934) used a Ba/Al alloy. 


; ! Tracings of oscillograms for certain Ba lines, for two different pulse lengths 
“Are given in Figures 3(a), (6). Figure 4(a)-(f) show reproductions to 3/5 full 


Current 
y 


co 


= Current = 
20 Ba 38924 53. Ba 4900A 
SZ Light. Emission 55 
ss 2 Shes 
se Y gs 
Be == Light Emission 
SS “iY we 
2 s bj WY .- y 
Xoo or 
= gee 0 
Oe 2) STAG ee Om2eI4 LG 
A Sec 
(a) (0) 


Figure 3. Oscillograms, light emission/time and current/time for sparks in air between 
Ba electrodes. 


}the current pulse and timing oscillation (2-8 psec. periodic time). 

| The light emitted by the sparks was shown to increase with time during the 
discharge and to reach its maximum intensity about 1 psec. after the current had 
begun to fall (Figure 2). The light is still intense when the current has fallen to 
izero, and diminishes gradually thereafter. 

The excitation temperature derived from the relative intensities of the 
above mentioned Mg triplets is plotted as a function of time in Figure 5. The 
results are discussed in §6 (from which, for these discharges, electron tempera- 
tures=excitation temperatures). 

It was found that ionized lines (Ball) reached their maximum strength of 
emission some 2-3 psec. later than neutral lines, from whatever element (Ba, Cd, 
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Figure 6. Comparison of Jight/time oscillograms for neutral and singly ionized lines. 


Mg) the latter were produced. The smal] peaks shown at the beginning of these 
and other oscillograms (Figures 4 and 6) are due to electrical pick-up. The effect 
was easily shown (since it arose from the current used to fire the 3-electrode 
triggering gap placed in series with the test gap) not to effect the oscillograms after 
the first microsecond. The spark currents are so high and the current pulses so 
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sharp and short in time that slight pick-up on the amplifiers proved very difficult 
to eliminate, since the overall amplification of the muliplier and amplifier system 
was 108-109. 

Spectra of Ba, Ca and Mg obtained with the short sparks are given in Figure 7 
(b), (c), (2) with Hg comparison spectra. The lines were identified from those of Hg 
by the use of standard interpolation methods. ‘The procedure for comparing line 
intensities was as follows for Figure 7 (a),Ca. The lowest exposure was taken first, 
and was obtained with an N.P.L. calibrated tungsten filament lamp (2700°K.) 
and with a Hilger step wedge mounted immediately before the spectrometer 
collimator slit. ‘The step wedge was removed and the Hg line 5461. super- 
imposed on the continuum; the other lines from the Hg lamp were absorbed by a 
suitable filter. "The next exposure was with the 4usec. Ca spark in air (35 psec. 
collimator slit), again with the step wedge, and was of 10 min. duration. In other 
cases the exposures ranged up to 70 min., the longest exposures being those taken 
with the 1 usec. discharges. The final exposure (1 min.) was with the Ca spark, 
without step wedge, and was taken to check the uniformity of illumination along 
the length of the slit. 

The distance of each microphotometered line was measured relative to the 
superimposed Hg line and corresponding points were then taken along the con- 
tinuum for density measurements. From the data thus obtained, using the six 
steps of the step filter for the meta] lines and the calibration continuum, the 
spectra] response of each plate and thus the relative line intensities were deter- 
mined. ‘The microphotometer technique followed normal practice. 

Since it appears (Hopwood and Craggs 1947) that the spectral response of the 
plates used in these experiments is different for the spark light (~1-10ysec. in 
duration), and for steady light, comparisons were only made between lines whose 
wavelengths differed by amounts (<400 a.) too small to cause appreciable errors. 
The existence of this phenomenon (intermittency effect on spectral response) is _ 
usually and justifiably ignored since it does not occur in most spectroscopic work. 


‘Table 1 

Bal (Kruithof 1943) 
Initial level (ev.) 3-80 3°79 3-16 3-16 3°15 3-14 
Wavelength (a.) 5778 5519 6063 5997 6019 6111 
Ag (relative) 6°55 3-41 1:95 0-84 1:00 3-06 

Ball (Mason 1938) 
Initial level (ev.) 2°5 5:22 5-22 5-66 
Wavelength (a.) 4934 4525 4900 3892 
Ag (relative) 1-00 0:95 1°55 3-41 

MglI (Kersten and Ornstein 1941) 
Initial level (ev.) 5-1 5-1 5-1 5-95 5:95 5:95 6:55 6°95, 
Wavelength (a.) 5184 5173 5167 3838 3832 3829 5528 4703 
Ag (relative) 50 29°7 9-9 o25 195 63-7 15 27°5 


Excitation temperatures with barium and magnesium electrodes were calcu- 
lated for equation (1) by substitution of the relative intensities determined and the 
appropriate Ag values, wavelengths and energy values. 

With magnesium the mean result was of the order of 10,000° k., but in the case 
of barium and in determinations on the barium lines carried out with barium— 
magnesium electrodes the fluctuations were very much greater. The excitation 
temperature determined from the values obtained for the ionized barium lines was. 
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with v4 correction) 1 :3 :5-2 respectively, while from intensity determinations 
using the above technique, the ratios were 1 : 3:2 :4 (again with v* correction). 

Further data showed that for the Cd triplet (4678, 4799, 5085 a.) the intensity 
ratios were 1 :2-4 :4:35 and 1 :2-7 :4-2 for two separate runs. For a Zn triplet 
(4680, 4722, 4811.) the ratios were 1 :2-6 :4:35. All these results include the 
vt correction. The theoretical values are 1 :3 :5 and Ornstein, Hengstum and 
Brinkman (1938) obtained for the above Cd triplet an experimental value of 
1 :3-4:4-2 with an arc discharge. It is generally considered that ratios having 
values so far from the 1 :3 : 5 value given by the Sum Rules indicate self-absorption 
for one or more of the lines. However, Mason (1938) found ratios of 1 : 3-32 : 5-37 
for the Zn triplet and considered this result to be so near the theoretical value that 
self-absorption could be ignored. Mason and others point out that the Sum 
Rules are not necessarily valid for these wide multiplet separations, so that it seems 
reasonable to suppose that self-absorption was not important in our experiments. 
‘This is borne out to a first approximation by the general consistency of the temper- 
atures given in table 2 for any particular set of conditions. Ornstein, Hengstum 
and Brinkman (1938) found ratios of 1 :3-7 :5-3 for the Cd 63S—5%P triplet and, 
it is interesting to note, ratios of 0-98 : 2-92 :5 for the multiplets 62Ds » ;-5 Ps 1,0. 
A close triplet (Zn) with components at 3018, 3036, 3072 a. gave ratios of 0-99 : 
3-13 :5 (Schuttevaer and Smit 1943) which is extremely close to the Sum Rule 
value (see also Smith 1945); if the available resolution had been higher in the 
present experiments, this triplet intensity ratio would have been measured. 
Further experiments with ultra-violet lines will be undertaken. The results given 
in Table 2 are probably not very accurate since some self-absorption may be 
present, on the basis of the above discussion. 


§5. EXPERIMENTAL RESULTS ON TRANSPORT OF METAL 
VAPOUR 


A brief investigation of some electrode effects is described in order to emphasize .- 
the difficulty and approximations in interpreting the results of Table 2. The 
results described in this section have a wider bearing since they should also be of 
interest to spectroscopists engaged in analytical work. 

Data on the transport of metal vapour in the discharge was obtained using two 
experimental methods : (a) observations on the spectrum of the electrode material 
in the discharge with an electron multiplier mounted in the focal plane of the 
spectrograph, and (b) photographing the discharge with a rotating mirror camera. 
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Figure 8. Collimating system for vapour transport measurements. 


In the first of these methods, the intensity-time variation of a number of lines 
in the spectrum of cadmium, which was the electrode material used, was obtained 
for different regions along the length of the spark channel. Cadmium was used for 
these multiplier experiments on metal transport since other metals ex 2n, Ba 
and Mg showed such large spark-to-spark variations that the oscillograms were 
useless. In order that such measurements could be carried out, a collimating 
system as shown in Figure 8 was used. G is the spark gap, and an image of the 
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Figure 11. Rotating mirror and still records for sparks in air between Ba ((a), (6), (c) 
9 usec. sparks, (d), (e), (f), (g) 4 usec. sparks), Mg ((h),(¢) 4 usec. sparks, (7), (k) 
9 usec. sparks), W ((L), (m) 9 usec. sparks), and Cd ((), (0) 4 usec. sparks) electrodes. 


ST a OT a, TE f 


$ ~ 7 = = 
. , 2 5 _ 
= a ~ ' 4 . 
5 : ao “ : & 4 
i ‘ 4 1x 5 ' 
: . . - i y * - > . 
- = a > = q P % oy 
. 2 i a 2 < m/ ; ’ 
ad > es i - E: _ : \ 
oar = [\p \ 4 
= =1 B ‘ 
», ‘ ' ¥ 
z . < . s 
. ’ or sa 
- 5 r= Ss " 5 ? 
* ; = “ ; : . j F ; 3 - ¥ " ; 
, 3 r 3 a i- 
; 2 - : a * = ‘ - 
> 7 ‘ : é a c oo -. 7 
“ . F “ae 
= :: . . ta 2 . \ ~ 
- . 4 - = ; . . 
- 
e ’ 
? 3 4 
a is 
: i } . = b . ‘ 
ve 7 a — Pre ap r 
* ‘ ‘e 
oat oe it < 
- a A , Le + , a @ 7 2 
’ 2 - 2. —p- = 
‘ 7 | > ~ ‘ > = 7 - 
< a = 7 - ~ - > ~ “ e 


Short Sparks in Air 57 


spark was focused by means of the lens L, into the plane of the screen S,. The 
screen carried a horizontal slit of variable width, and was mounted on a catheto- 
meter, so that easy vertical movement and measurement of the slit position in the 
image was possible. Thus light from a small region of the image of the spark was 
allowed to pass through the collimating system into the spectrograph slit. An 
image ofthe slit S,; was focused by means of the lens Ly, which was mounted directly 
before the spectrograph slit, on to the collimator lens. This was to ensure that all 
light which passed through the slit passed through the optical system of the 
spectrograph, and no false background was produced by reflection from the 
collimator walls. 

In the experiments performed the length of the spark was about 6 mm. and was 
slightly magnified by the lens L, so that an image ot length 8 mm. was produced. 
‘The width of the slit S was 0-1 mm., and observations of the light-time variation 
were made at intervals of 1 mm. along the length of the image. 

Measurements were carried out on the two Cd lines 5085. and 5379 a. for 
currents of about 120 amp. (duration 4 usec.), and about 30 amp. (duration 10 psec.). 
‘The observations were repeated a number of times over a period of weeks. 

Figure 9 (a), (b) show typical oscillograms taken for a 10 usec. spark at 1mm. 
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Figure 9. Light/time oscillograms with cadmium electrode sparks. 
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Figure 10. Transport velocity of cadmium vapour in spark gaps. 


intervals along the spark image. The first part of the oscillograms, being slightly 
affected by pick-up (Figure 4(d)), has been omitted. It can be seen that the onset 
of radiation takes place progressively later, and the radiation becomes weaker as 
the distance from the cathode increases. In some of the oscillograms it was 
observed that after having moved some distance along the channel the radiation 
became stronger again and occurred earlier in time, suggesting that some metal 
vapour was being given off at the anode. This was confirmed by the rotating 
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mirror photographs; with Cd electrodes this anode radiation is much weaker than 
that originating at the cathode. The spark to spark variation in the intensity of 
the light originating from a small region of the discharge is much greater than that 
in the integrated light at any one wavelength. 

It is possible from’ these results to determine the velocity of transport of the 
metal vapour through the gap (Figure 10 (a) and(d)), and it appears from the curve 
that the velocity is approximately 1-1 x 10°cm/sec. for the 4usec. discharge, and 
that in the early parts of the 10 psec. discharge (first 4 usec.) the velocity is 
approximately the same, but falls steadily to about 4 x 104cm/sec. after about 
8 psec. 

Typical rotating mirror and still records for single sparks, using the same 
current waves as before, are shown for various arrangements in Figure 11. In all 
the photographs the anode is the upper electrode; the peak current in each case 
was 110 amp, the still records are Figure 11 (c), (g), (2), (R), (m) and (0). 

There are many inexplicable features in the records of Figure 11. The 
extraordinary variation in the amount of material volatilized for each spark with Ba, 
which accounts for the great fluctuations in size of the electron multiplier light 
oscillograms, is to be compared with the steadiness of a discharge with Mg or Cd 
electrodes. The light emission from metal vapour emitted from either electrodes 
is approximately equal with Ba, Mg or W electrodes, but with Cd most of the 
vapour originates from the cathode. In the latter case the vapour speed is much 
higher (as measured with the multiplier) than with Ba or Mg photography, and it 
seems likely that several mechanisms of vapour emission are operative (§6). 
With Ba for example, there are several clouds of vapour which appear to move very 
rapidly into the mid-gap region and then lose their luminosity without further 
movement (f). There are slowly moving vapour clouds noticeable with Ba (e) 
and also with Mg (7), but the latter does not seem to show rapidly moving clouds. 
The afterglow in the channel (due possibly to nitrogen) is seen in (/), W electrodes, _ 
since the width of an infinitely thin 4 usec. discharge should be only ~1 mm. with 
the speed of mirror rotation used. The sharp channel can be seen on the left of 
the spark image in (h) and (I). 

Further work on the mechanism of metal vapour emission in air and other 
gases is planned, using long and short discharges, but the present relevance of the 
records (Figures 10 and 11) is to show that the excitation temperatures as found 
for these discharges by the methods described are only time and space means. 

Vapour speeds deduced from such records as those in Figure 11 for various 
electrode materials and current pulse lengths are as follows: Mg, 4ysec., 
2-4. x 104 cm/sec. from the cathode, 2:6 x 104 cm/sec. from the anode; Mg, 
9 psec., 2-6 x 10% cm/sec. from the cathode, 1-2 x 104 cm/sec. from the anode; Cd, 
4 psec., 2-2 x 104 cm/sec. from the cathode; Ba, 9 psec., 1-3 x 104 cm/sec. from 
the cathode. 


§6. DISCUSSION OF RESULTS 


The excitation temperatures (‘Table 2) show that a statistical electron distri- 
bution, giving the same “temperature”? for all levels, does not exist. Langstroth 
and McRae (1938 a, b) arrived at the same conclusion for fairly short sparks. The 
meaning of the temperatures of Table2 needs clarification and reference should be 
made to the work of Ladenburg (1933) who points out that in cases where a thermal 
equilibrium. does not exist a statistical equilibrium between excited atoms and 
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electrons may befound. The electron velocity distribution should be Maxwellian 
and the excitation to and destruction of excited states must be due largely to: 
electron impact.* 

Ladenburg shows that the Klein-Rosseland (1921) relation 

Bi Ae \aen EOL) 3 ee ED (2) 

may be applied to plasma conditions without the necessity for assuming a statistical 
equilibrium between electrons and excited atoms, since the S’s and g’s have values. 
independent of an equilibrium; g, and g, are the statistical weights of the two 
energy levels, FE; and E,. £” and E’ are the energies of the exciting electrons 
(before and after collisions) so that E” >E,—E, and E’ = E” —(E,,— E;)..  Syj(E’) 
is the probability of destruction of an excited atom ina collision of the second kind’ 
with an electron. Ifthe current density in the plasma is very high (as it is in spark. 
channels), so that many such collisions occur, then the number of excited atoms 
destroyed per second is equal to the number produced, other processes being 
considered negligible; Ladenburg then finally finds that 


NyINj=(eele)exp(—(Ex—E)/RTI, 0 sees (3) 
i.e. the ratio of populations of the two atomic levels corresponds to a statistical. 
equilibrium at the electron temperature. ‘The temperature of the normal atoms 
can be much lower than 7,. An excellent review of these and other matters was 
published by Llewellyn Jones (1945). 

It appears from the results (Table 2) that the same 7, (about 10,000°) applies to 
the MglI levels studied, independently of current pulse length. For Bal,. 
greater variations in T, for the different levels are found and 7, (Bal) is greater 
for sparks between Ba electrodes than with Ba/Mg electrodes. T, for the Ball 
levels is consistently higher than for Bal, reaching 20,000°.. The relation between 
the nature of the electrodes and T, is probably complicated and might be a function 
of electrode melting point, vapour pressure, and excitation level energies etc. 

The persistence of luminosity of the metal vapour after the current had fallen 
to zero (Figures5 and 11) is sufficient to prevent any information being obtained on 
the rate of loss of electronic energy in a currentless plasma at high pressure.. 
The effect may be merely due to the finite lives of the excited states in question, 
but might be due to yet another complication, i.e. the excitation of metal atoms by 
collisions of the second kind with metastable atoms produced in the air during the 
earlier part of the discharge. The work of Manley and Duffendack (1935) and_ 
Duffendack and Thomson (1933) on line enhancement, e.g. of magnesium by neon, 
demonstrates the feasibility of this explanation. Bay and Steiner (1930) found, 
for example, that the metastable (a2) nitrogen molecule could excite a mercury 
level (23S) at 7-7 volts. Further work with different gases will, it is hoped, give 

* One of the referees of this paper has kindly brought to our notice the review article by Rompe- 
and Steenbeck (1939) who discuss in detail “non-isothermal” plasmas. They point out (as 
have many other authors) that in many discharges, the electron, excitation and gas temperatures 
may be widely different although at higher pressures (101° atoms/cm* or more) the electron temper- 
ature is only slightly less than the gas temperatures (say 20-100° less at 6000°) in steady maintained 
discharges (arcs). ‘Rompe and Steenbeck do not treat transient conditions (cf. the work of Mann- 
kopff and Paetz 1947). It seems reasonable to suggest, from the work of Ladenburg (1 933) and of 
Mannkopff and Paetz that, in the sparks discussed in the present paper, thermal equilibrium will 
not be attained but that the excitation and electron temperatures will probably be almost equal 
(relative to a particular excited state) and that the several microseconds available during current flow - 
should be enough for statistical equilibrium to be reached amongst the electrons. It seems that 
the excitation temperatures measured relative to the various energy levels. are different (Table 2): 


and further work is planned to investigate this and other matters in more detail. Experiments to- 
measure electron temperatures separately are also in progress. 
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further information on this point, but the advantages of the new pulse techniques 
should enable the experiments to be simplified, and render them capable of adding 
new data to the earlier results of Duffendack and his collaborators. 

The mean life (r,,) of an excited state 7 is derived from a knowledge of Einstein 
A’s since 7,,=1/XA,» where m is a lower state. 

The A’s must be absolute and not relative or as those given in Table 1, and such 
-data seem to be lacking. | Other methods are available for finding 7 (Mitchell and 
Zemansky 1934, pp. 145 et seq.), but the 7 values given are not for those states used 
in the present work. It is of interest to note that 7(Zn, 3076.) is ~10 psec. and 
7(Cd, 3261 a.) is ~2 pusec., so that true lives of the order of those shown for example 
in Figure 4 are not unreasonable. 

It is assumed that in the time of current duration (1 sec. or more) there have 
been sufficient collisions to enable the equilibrium between electron and excited 
atoms to be established. This seems reasonable in view of considerations 
involving the number of electronic collisions in 1 psec. in spark channel conditions 
(Ladenburg 1933, Llewellyn Jones 1945, Margenau et al. 1946). ‘This number of 
collisions is probably about 10°, and the time required for the electrons in the 
channel to reach an equilibrium energy distribution should therefore be not more 
than about 1 wsec. The mean life of the excited states must be then much less 
than 1 psec. for the Ladenburg equilibrium to be established, in which case the 
long afterglows could not be explained (see above) by the persistence of the excited 
atoms produced during the time of current flow.* 

The mechanism of metal vapour emission at the electrodes may affect the values 
of T,. Langstroth and McRae (1938b) consider that diffusion or transport as 
charged particles could not account for the high vapour speeds (~10° cm/sec.) and 
suggest that an explosive process at the electrodes might be operative. Haynes 
(1946) in a short abstract has described in outline some work on the velocity of 
‘mercury vapour jets in hydrogen sparks of microsecond duration, using electron 
multipliers. Haynes concludes, without detailed analysis, that the high speed 
clouds or jets might be caused by the positive ions and electrons acquiring momen- 
tum in the anode and cathode drops respectively. This is an obvious general 
interpretation but it seems unlikely to hold for al! the cases shown in Figure 11. 
"The difference in behaviour for cadmium and magnesium electrodes is difficult to 
explain on such grounds, since (Mg, Figure 11(j)) it is sometimes found that the 
anode and cathode clouds are very similar in appearance and speed of propagation. 
Llewellyn Jones (1946) has produced strong evidence and arguments in favour of a 
simple evaporation process, which is supported by independent work (Craggs, 
Haine and Meek 1947). In this case it is still difficult to account for the different 
behaviour of the various metals. 

The movement of a cloud from the anode might be due to positive ions. 
‘Tyndall (1938, pp. 67 et seg.) shows that the drift velocity W of a positive ion is not 
proportional to X/p (field-strength/pressure) for values of X/p>10v/cm/mm. Hg. 
The mobilities for Mgt, Cdt+t, and Bat in nitrogen are about 3, 2-3, and 
2:2 cm/sec/v/cm. respectively at N.T.P. and the variation with temperatures up to 
500°k. is negligible at constant density (Tyndall 1938, pp. 45, 58 et seg.). The 


* Mannkopff and Paetz (1947) have studied the relaxation time for the electron temperature in 
an ionized gas, and for atmospheric conditions find this time, i.e. the time for a Boltzmann equili- 
brium to be established, to be much less than 1 usec.; they also state that at such high pressures the 
difference between electron and excitation temperature is negligible. More experimental data 


‘seem to be required and work to this end has been in progress at the University of Liverpool for 
some time. ‘ 
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initial gradient in the present spark gap is approximately 40 kv/cm. so that the ions. 
mentioned could have drift velocities of approximately 8 x 105cm/sec. which are 
adequate to account qualitatively for the speed of the anode vapour clouds. 
Further speculations, particularly on the cathode clouds, seem unjustified. The 
data on massive negative ions is extremely scanty. 


§7. CONCLUSIONS 

Excitation temperatures in 110-amp. sparks lasting 1-8 psec. have been 
measured with various electrode metals and values varying between 3000°k. 
and 20,000°k. have been found. There is no clearly defined equilibrium between. 
the electrons and excited atoms of all levels in the spark channels. 

The results for these discharges, and those of the few earlier workers experi-- 
menting with similar discharges, must be interpreted with caution since the 
mechanisms by which the luminous vapour clouds are formed at the electrodes do. 
not appear to be simple. . 
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§1. INTRODUCTION : 


“YN a previous paper (van Itterbeek and van Doninck 1946) we published a 
series of measurements on the velocity of sound in gas mixtures at low temper- 

atures, the components being monatomic gases or hydrogen which possess 
“no rotational energy at low temperatures. 

In that paper we developed a theoretical method to compute the interaction 
term between the components, starting from measurements on the dependency of 
‘the velocity of sound as a function of pressure at fixed temperatures. 

In the present paper we report on measurements of the same kind, but with 
components for which the rotational energy must be considered. So measure- 
‘ments are made on H,—-N,, H,-O,, H,-CO, He-N,, He—-O, between 90°. and 
65° K. 

Just as for our previous measurements we find that the second virial coefficient 
B for the mixtures can be represented by means of the equation 


B=B,x? + Bx(1 —x) + B,(1—x), 


_B, and B, being the virial coefficients for the components. The mixtures H,—N, 
-are however an exception to this rule; in this case good agreement could only be 
-obtained with the measurements by means of the following expression for B: 


B= By,x? + Bx3(1 —x) + By (1 —x)?. 


We suppose that the reason for this deviation is the great chemical affinity 
between the hydrogen and nitrogen molecules as it appears from the formation of 
-ammonia. ‘This fact led us to measure also the mixtures H,-O, and H,-CO. 
But as for the former mixtures, no deviations appear from equation (1). 

The experimental arrangement has been described already, in part in the paper 
already referred to and in another paper in which this technique was first used 
(van Itterbeek and van Paemel 1938). It consists in essence of a quartz oscillator 
‘of frequency 523-78 kc/s. emitting ultrasonic radiation into a gas contained in a 
vessel in which the pressure can be varied. The movement of a reflector mounted 
on a micrometer screw enabled wavelengths in a gas to be measured by the reaction 
of the returning radiation falling on the emitter. 
ments to one thousandth of a millimetre. The frequency of the quartz was 
measured at the boiling point of oxygen. A further measurement in liquid 
hydrogen indicated that variation of frequency of the quartz with temperature was 


not so important as possible changes in the temperature of the gas during a series, 
which constituted the greatest source of inaccuracy. 


The screw permitted measure- 
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§2. THEORETICAL FORMULAE 


The velocity of sound as a function of pressure in the mixture is given by the 
equation ; 


W=W,(1+pS/RT), 


where W, is the velocity of sound in the gas mixture in the ideal state (p =0), 
which can easily be computed from the concentrations and the specific heats of 
the components, R is the gas constant, T the absolute temperature, and 
TdB IRE ake; 
S=B+ —-— + ——— —_,—=8a 
aT” ALD aT? By 

A is a numerical factor which can be computed from the specific heats and the 
concentration. Sis computed by means of the method of least squares from the 
experimental curves represented analytically by equation (3). 

For the case of the H,-N,. mixtures we found a good agreement with the 
experimental data by using the equation 


S=S'yxe+ax(1—~x)+S’y(l—-xF ...... (5) 
and 
; T dBy, IR eee . 
ear 20, e) ate ©) 
and an analoguus expression for S’y,. 
Further as 
1 dp ie d?B 
ce 7 
seg oes de (7) 
and Nn =A +AQ(1 — x); Ap= (CH), <o/My. 
For all the other mixtures we found 
Sy Mean wel x) aaa es (8) 
As for our previous measurements we found that the expressions 
a=fS —S,'x'—S,/1—x)]/e8(l—x)  — ...... (9) 
-and 
GS = )S,x2—S,(1—x)")fa(l—2)' © ees (10) 
are practically independent of concentration, which proves that the expression 
T dp 1 SoG 


as ees ae ap SSS 

AX dT 2XA+1) dT? 
is small. Thus with good approximation we can write «=f and obtain some 
indication of the interaction between the components of the mixture. 


(i) Measurements on Hydrogen Mixtures 


(a) H,-N,. ‘The velocity of sound as a tunction of pressure between 0-1 and 
0-8 atmos. at liquid oxygen temperatures (90° K., 85°K., 80° K. and 75°.) and for 
different concentrations has been measured. From these measurements we 
have computed the values of W, and S (see Figure 1, which gives typical data at 
80°x.,and Table1). In this table we have compared the experimental values for 
W, with the values obtained theoretically for the mixture of the perfect gases. 
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In Table 2 are indicated the values for «=f computed by using equation (9). 
It appears that « is practically independent of concentration (cf. Figure 2). 
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2 en ee 


-10° s 


i96— 


Lee P) 90-9%N,-9-1%,H. 
Se ed 


1g9—— 
183 


178 —L 
0 0:2 0-4 0-6 0:8 100 80 60 40 20 
p (atm) | #N, 
Figure 1. Figure 2. 


(b) H,-O,. Similar experiments have been made on H,—O, mixtures; the 
results of the calculations are indicated in Table 3 and Figure 3, and the values of 
a orf are computed by means of equation (10) and given in Table 4 and Figure 4. 


———— 


H,-0, 7=80'k 


27-9L,0,- 72:1, Hy 
50-06% 0,-49-947%H, 
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Figure 3. Figure 4, 
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Velocity of Sound at Low Temperature 


Table 1 
Concentration 
YN,  %He 
100 0) 

91-3 8-7 
74-6 25-4 
42°8 57:2 
21:2 78°8 
ee) 100 
100 0 
88-2 11-8 
76-4 23-6 
62:0 38-0 
43-8 56:2 
16:0 84-0 
0 100 
100 0) 
90-9 9-1 
84-7 15:3 
60:4 39°6 
45°5 54:5 
24:8 75-2 
@) 100 
100 0 
90-9 9-1 
75-9 24-1 
62-0 38-0 
44-1 55-9 
34-1 65-9 
0 100 

Table 3 
Concentration 
%O2 WH, 
100 0 

81-2 18:8 
71°5 28°5 
50:8 49-2 
28-4 71-6 
0 100 
100 @) 
83-0 17-0 
64:7 35:3 
50-0 50:0 
33-3 66°7 
0 100 
100 0) 
84-0 16:0 
63:1 36°9 
50:1 49-9 
27-9 72:1 
0 100 


W, (exp.) 
(m/sec.) 
193-8 
205°5 
224°5 
291-7 
392-1 
- 773-0 
188-3 

200-4 - 
2157 
238-4 
280-4 
422-2 
PSS 
183-1 
192-2 
199-0 
235-2 
268 :9 
348-4 
731-4 


176°3 
186-0 
204-2 
224-0 
263°3 
294:7 
707-0 


W, (exp.) 
(m/sec.) 
181°0 
200:2 
214-7 
253:2 
329-1 
773-0 
175°6 
193-4 
219-0 
248-2 
299-2 
Wales 
171:0 
186-6 
215:2 
240-4 
310-4 
731°4 
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H.-N, Mixtures Table 2 
: W, (theor.) oe Ag , 
NOES SS (aiisee) INCRE) AN Oe s< 
—4-81 193-8 85 90 5) 
— 4-03 205-4 80 5:9 
—2:57 224°3 70 5:8 
—0-93 292-0 60 6-0 
—0-94 392-4 50 6-4 
0:29 40 6:4 
—5-96 188-1 80 90 8-1 
—4-31 200-2 80 8-0 
—2-94 215-8 70 7:6 
—2:16 238-2 60 8-0 
—1-:17 280-4 50 8-2 
23 422-1 40 oe 
0:24 750.90 9-0 
—6:68 182°8 80 &-9 
—5-03 191-9 70 9-2 
—4-28 198-7 60 9-8 
—9-64 234-9 50 10-0 
—1:57 268-7 40 10:0 
—0-76 348-2 
0-17 
—8:21 176:6 
—6:32 LSDe7, 
—4-24 204-3 
—3:37 224-1 
—1-91 263-0 
—1-43 294-7 
0-14 
H,—O, Mixtures Table 4 
W, (theor.) 3 x Be 
10?x S$ Gifseed fa (GeKe) Or Mt xe 
— 5-90 90 80 —2-44 
— 4:31 200°1 70 —2:83 
— 3.66 214-4 60 —2:-98 
— 2:26. ASI 50 —3-:05 
— 0:87 328°3 490 —3-06 
0:29 30 —3-06 
= VION 85 80 —4-78 
— 6:07 193°5 70 —4-82 
— 4:59 219-1 60 —4-83 
— 3-10 248-1 50 —4-83 
— 1:97 299-0 40 —4-83 
0:24 30 —4-82 
—10-60 80 80 —5-2 
— 8:95 186-4 70 —6:0 
— 5:98 215.3 60 —64 
— 4-72 240°5 50 —7-2 
— 2-40 310°5 40 —7:0 
0-17 30 —7-2 


6:0 


8:0: 


9-6 


10° x a 
average 


= P39 


4.32 
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(c) H,-CO. Figure 5 represents typical experimental data for the velocity ; 
Figure 6 and Table 5 give the computed values for S. 


H,-CO | 7'=80'k. 


22:74.00-77-30H, 


6, » in) 

> 37 64,00 -62-4/%H, a 
3 3 ' 
290 

= 

7k 

Ss 
~~ 


62:6%, C0 - 37-4/,H, 


100 80 60 40 ry 
460 
Figure 5. Figure 6. 
Table 5. H,-CO Measurements 
Bee Concentration W, (exp.) 3 W, (theor.) 
CES) %CO oH, GES to (apece) 
90 100 0 L937 
Tie? 28:8 229-6 —3-32 229-6 
46°3 S307) 282-1 —1-74 281-9 
29:7 70:3 342-7 —0-98 342°5 
7S 82-7 423-3 —0°55 423-0 
85 100 0 188-2 —5-92 188-0 
88-3 ilalay/ 200-2 —5-10 200-4 
53:8 46:2 254-7 —2-59 255-0 
Bit 62-9 302:3 —1-70 302-3 
DUES a7 S185 —0-84 Set 
80 100 0 183-0 
80-0 20-0 204-1 —4-75 204-3 
62:6 37-4 230°6 —3°:51 230-4 
37:6 62-4 291-9 —1-83 291°8 
227, US 360°3 —1-08 360-0 
75 100 0 7D: 
87-4 12-6 189-4 —6°37 189-3 
68-6 31-4 PUGS —=5-00 2133 
42-3 S707 268°1 —2-82 268-0 
We 82-3 382-9 =1-07 383-0 


For these measurements « and 8 could not be computed because the virial 


coefficients of CO are not known. We found however that JS as a function of the 
concentration could be expressed by means of the equation 


is) —— Soox? + ax(1 =e) + Sy — x), 


Soo was taken from our earlier measurements (van Itterbeek and van Doninck 
1943). The values of « can however be considered as approximate values for B. 
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Phese are 10x = ~3-29, 3-96, —426, —5-43 for T=90, 85, 80, 75°x. respec- 


SvVeL¥. 


(33) Measurements on Helium Mixtures 


© He-N, Figure 7 gives W as a function of p and the concentration at 
wx In Table6 oe of S and in Table 7 the values of z or 8 


‘cf. Figure 8). 


asic 


; 


LIM CMe 
- ; E2479, 4T Gite 
/ 
) ) 2 
———__| ee = 
& s BEM, D4i He 
= 4 ; 3 ss 25 
1 : ; mass 
; | i / ; 
| NIM, BT Me / 
; e O- | 
i | 
£2 ge te 62 190 9 
p Gi) ALP 
Figure 7. Figure 8. 
Table 6 He-N, Mixtures Table 7 
Te Concentration Wg exp.) W, (theor.) 5 x eNO? 500 
EOD oN, He (mice) “*S misec) 7%) Gn, 17% average 
WD B+ 96 034 —413 203-5 9 80 —10 
$59 1414 2104. -—351 2106 70 —10 
627 373 WOO -—212 1407 oO +13 4, 
428 372 B19 -095 282-4 30 ~—1+4 
9 1090 5595 0-50 559-2 49 —13 
3 873 127 WH  -490 2008 Oe Ms 
617 383 2357 —271 235-4 30 80 —26 
«424A 3576 SA 1-43 275-5 70 —2-4 
296 W4 3154 —-069 315-2 Ne ee 
9 02 54249 0-38 3427 ; 30 —2-9 
0 NZ 87 1908  -563 1908 5 iene eed 
BG 4 M45 443 047 SB) ee 
324. 476 «=—52 250 245-1 75 80 —40 
339 = 65-1 288-2 1-37 288-2 Di 37 
9 2 52649 0-28 527-1 60 —3-9 
eae e 
i 2 7-4 1835 —719 183-4 
BF 233 M7 -525 m29. a 
6H 324 BOB -270 270 aa 
322 618 WoT 167 2865 
9 0 359949 0-19 3102 
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(0) He-O,. Figure 9 represents typical experimental results of W as a | 


function of pressure. 


Table 8 contains the values of S (Figure 10) computed — 
from the measurements, and Table 9 values of «. | 


Table 8 
o.., Concentration W, (exp.) 
NS %O, %He (m/sec.) 
90 100 0 181°0 
89-1 10-9 191:8 
64-3 ye7/ DP aiche! 
41-7 58:3 269-8 
22-6 77:4. 338-3 
85 100 0 175°6 
89-1 10:9 186-0 
64-3 35°77 2074 
41-7 58-3 262:1 
22-6 77-4 328-9 
80 100 0 171:0 
89-1 10-9 180-4 
64:3 SI0)97/ 210-6 
41-7 58:3. 254-3 
22-6 Tt BAS 7 
= =) 
| “He-0, 7'=80'k. 
319 


Sit. 


22-640,-77-44He 


255 | 
© 41:77/0,-58:34He 
3 253 — 


89-17,0,-10:9% He 


04 


P (atm) 
Figure 9. 


0-6 


He-O, Mixtures 


10°x S 


W, (theor). 
(m/sec.) 


191-6 
22326 
269-7 
338-4 


186-2 
PATA 
262-1 
326°7 


180-6 
210:8 
254-1 
319-0 


Table 9 
° x 3 
De) Oz 10? x a 
90 80  —1:03 
70  —1:33 
60 —1:50 
50 —1:56 
40 —1:58 
30 —2-00 
85 80 —1-99 
70 —2-10 
60 —2:-29 
50 —2:-24 
40 —2-42 
30 —2:-52 
80 80 —3-32 
70 —3-42 
60 —2-89 
50 —3:-24 
40 —3-65 
30 —3-52 


103 x @ 
average 


—15 


—2:26 


3.44 


§3. DISCUSSION 


fo 92 
Figure 10. 


As we have said in the introduction, from our experimental results we can 
conclude that the second virial coefficient for the different mixtures can be repre- 


sented by means of an equation of the type (1). 


The mixtures H.-N, make an 
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exception to this rule; for them other powers than the second must be chosen for 
the concentration. ae - 
In Table 10 we have summarized the different values of « together with the 
second virial coefficients for the components. 


Table 10 
Gas mixtures (AR) 10°xa OZ eB 10° x By 
7 H.-N. 85 6:0 —0:3 = 9-3 
80 8:0 —0-4 —10:3 
75 9-6 —0°6 —11:6, 
H,-O, 90 —2°-9 —0:2 = 9-3 
. 85 —4-8 —0°3 — 10-0. 
os 80 —6°5 —0-4 — 11-0 
H,-CO 90 —3°3 
85 —3-9 
80 —4-6 
75 —6-4 
He-N, 90 , —1:2 0:5 — 84 
80 —2°-7 0:5 — 10:3 
75 —4-0 0:5 —11-6 
He-O, 90 —1°5 0:5 — 93 
85 — 2:26 0:5 — 10:0 
80. — 3-34 0:5 —11-0 


From this table we see that there is a complete difference between the values 
obtained for the H.-N, mixtures and the other ones. 
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Ultrasonic Dispersion in Organic Liquids 


By A. SCHALLAMACH 
The British Rubber Producers’ Research Association, Welwyn Garden City, Herts 


MS. received 3rd May 1948 


ABSTRACT. The sound velocities v have been measured, at frequencies f between 
about 700 and 3,000 kc/s., and at temperatures between about +15° and —80° c., of the 
following liquids : geraniol, di-dihydrocitronellyl ether, castor oil and iso-butyl alcohol. 
There is evidence of a positive dispersion (@v/@f>0) in the first three liquids at low 
temperatures. A negative dispersion has been observed in geraniol and di-dihydrocitron- 
ellyl ether at higher temperatures, and in iso-butyl alcohol over the whole temperature 
range. ; 

The positive dispersions are discussed on the basis of the classical hydrodynamic theory, 
and the implications of a complex viscosity are examined. 


$1. INTRODUCTION 
N classical hydrodynamics, the amplitude of a plane sound wave in a viscous _ 
] liquid is given by (Lamb 1945) 


U=U, exp (Jott mx), 9” eee (1) 


m= —jw//(k/p))(1+4jan/3x)y*, sa. ee (2) 
where « is the direction of propagation, and w the pulsatance of the wave in a liquid 
of density p, viscosity 7, and adiabatic bulk modulus x. The sound velocity is © 
found from the relation 


with 


hos —I(m)\i@s>—— . ok aueeheee (3) 
If wy/k<1, one has, on expanding equation (2) 


wo [ace 2 (ae a 


where 
Opa lily. a) ae seen ease (5) 


Equation (4) shows that the classical theory predicts a dispersion of the sound 
velocity if the product of frequency and viscosity becomes comparable with the 
bulk modulus. Ordinarily, this effect is negligible and the standard text-books 
give therefore equation (5) for the sound velocity. 

Apart from the viscosity effect, no dispersion should be found as long as the 
mechanical properties of the liquid do not depend on the frequency. The bulk of 
the experimental evidence indicates that there is no ultrasonic dispersion in liquids 
up to about 100 Mc/s. (Parthasarathy 1936, Bar 1938, Krishnan 1939, 1941, 
Bhagavantam and Rao 1946). Exceptions are castor oil for which Zachoval (1939) 
reported a rise of the sound velocity with increasing frequency, and carbon 
disulphide where, according to Te-T’chao (1946), the velocity drops with increasing 
frequency. 

The sound velocity measurements at low temperatures reported in this paper 
have been carried out in certain long-chain liquids which do not crystallize at these 
temperatures but eventually set to a glass. It was thought that the continuous 


° 
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transition to a glass might affect the ultrasonic properties of these liquids 
even at only moderately high frequencies. At the same time, regions of 
temperature and frequency where dielectric relaxation occurs would be covered 
in these experiments, so that it would be possible to compare these two effects 
and to see if there were any connection between them. 


§2. EXPERIMENTAL 
Two suitable liquids used in this investigation are geraniol 
(CH) GHe (CH),..C(CH,)).CH.CH,..0H 
and di-dihydrocitronellyl ether 
(Gree Cl. ChH.).: CH(CH.)t (CH). 0. 
The relevant data for their purity were given in an earlier publication (Schallamach 
1946). Since a preliminary note on this work appeared (Schallamach 1948) we 
experimented, in addition to these long chain liquids, with castor oil of medicinal 


quality “cold drawn’’. For comparison, the sound velocities in iso-butyl 
alcohol were determined. The sample used was a commercial product. 


Inches 


Figure 1. The interferometer. 1, crystal; 
2, ebonite holder; 3, spring; 4, high 
tension lead; 5, reflector. 


The measurements were carried out in a sonic interferometer specially 
adapted for low temperatures. The essential part of this instrument is shown in 
Figure 1. The X-cut circular quartz with silvered electrodes was pressed against 
the rim of an orifice in the bottom of the liquid container, without an intervening 
membrane. In this way, good sound transmission was assured at all temperatures. 
In the course of two days, ample time for an experiment, a few drops oozed out of 
the vessel but this did not affect the measurements. As no packing material 
suggested itself for the other joints, metal-to-metal contacts were relied upon for 
tightness: the bold rings on one set of flanges pressed against a flat part of the other 
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flanges served this purpose well. The reflector was adjusted by a micrometer 
arrangement about 10 inches above the interferometer vessel which was immersed 
in a well stirred acetone-Drikold bath and maintained at an internal nitrogen 
pressure of about 10 inches of water. The temperature was measured by a 
thermocouple in contact with the interferometer and was accurate to +0-1°c. 

The quartz crystal was wired in parallel with a resonant circuit which was 
loosely coupled to an oscillator, and the voltage across the crystal was measured by a 
valve voltmeter. It is easily seen that this voltage indicates the equivalent 
impedance of the system crystal-liquid column. Hubbard (1931) has shown that 
this impedance consists of two series resonant circuits, connected again in series. 
One of the resonant circuits represents the crystal, and the other, representing the 
liquid column, is found to be formally identical with the input impedance of a 
transmission line. As long as the damping of the sound waves in the liquid is not 
too great the impedance has a minimum whenever the position of the reflector 
coincides with a node of the particle velocity. As the reflector is moved through 
the liquid sharp decreases in voltage across the crystal are observed at positions of 
the reflector half a wavelength apart. 

This simple relation between wavelength and minimum of the impedance 1s no 
longer strictly valid when the damping becomes strong. It is then necessary to 
evaluate the exact condition for impedance minima which is 


tanh (2al) = —(2m/Ax) tan (4771/A), 


where /is the distance between the crystal surface and thereflector, Athe wavelength 
of the sound and « the amplitude absorption coefficient. Abadie (1946) has given 
graphical solutions of the above equation. The distances between impedance 
minima are now less than half a wavelength, but this effect is serious only for the 
first half-wave and quickly loses importance as the distance between quartz and 
reflector increases. Actually, the construction of our interferometer did not 
allow measurements of the first few half-waves to be made and the error is small for 
this reason alone. ‘The greatest error likely to arise from this source has been 
estimated for di-dihydrocitronellyl ether at —80°c. The viscosity is known 
(Schallamach 1946), and the absorption coefficient « was calculated by means of the 
classical equation « = R(m) where m is given by equation (2). It was found that at 
2,102 kc/s. the error is less than 0-01%%, and that at 689 kc/s. less than0-1%. This 
value is, at —80°c., small compared with the error due to other sources and was 
therefore neglected. 

The accuracy of the wavelength determination depended mainly on the distance 
over which a standing wave could be established with certainty, and it became 
progressively worse as the temperature fell and the attenuation increased. Near 
room temperature where the damping is unimportant the wavelengths were 
accurate to +0-05%, and this accuracy was maintained in the long chain liquids 
down to about — 60° c. where the error began to rise until it reached +0-35% at 
about —80°c. 

The crystals were excited in their fundamental frequency and its third har- 
monic. ‘The highest frequency used was approximately 3,000kc/s. but most 
measurements were carried out at about 700 and 2,000kc/s. Using a heterodyne 
wavemeter the frequency could be measured to an accuracy within +0:05 ise 


‘Thus, the sound velocity determined as a product of frequency and wavelength was 
accurate to between +0:1% to +0-4%. 
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As a check on the apparatus, the sound velocity in water was determined at 
two frequencies. Water was chosen for this purpose because it represents the best 
established case of a liquid with no dispersion. ‘The values obtained (in m/sec.) 
were: 


‘Temperature ° c. 23 16-6 
Frequency | 689 1489 1468 

(kc/s.) 2102 1491 1470 
Willard (1947) 1493 1476 


‘This table shows that there is agreement between the measurements at different 
frequencies, and reasonable agreement with the values obtained from an empirical 
equation given by Willard. 


§3; RESULTS 


Figures 2 and 3 give the results of the last of a number of experiments with 
geraniol and di-dihydrocitronellyl ether; improvement in technique during these 
experiments did not materially affect the results. As will be seen, both graphs 
show the same type of curves: at any one frequency the sound velocity increases 
with decreasing temperature. ‘This increase is almost linear at higher temper- 
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Figure 3. Sound velocity in di-dihydro- 
citronellyl ether, vs temperature. 


© 689 ke/s. ; x 2,102 ke/s. 


Figure 2. Sound velocity in geraniol, 
vs temperature. © 689 kc/s.; 

x 2,102 kc/s.; A 2,943 kc/s. 
atures, and in the same region the sound velocity has a small negative frequency co- 
efficient. For example,comparing measurements in geraniol at 689 and 2,102kc/s., 
the relative change at 10° c. is 0:2%, which is near the experimental error, but 
at —50°c. it increases to 0-65°%%. From the number of points which give a 
consistent result, this effect appears to be real and will be referred to as negative 
dispersion. At lower temperatures, the temperature coefficient becomes con- 
siderably larger and the curves cross each other so that the frequency coefficient 
changes its sign. This positive dispersion increases as the temperature falls and 
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at —80°c. the relative difference between 689 and 2,102 kc/s. is about 3% in both. 
liquids. We should mention here that extending the temperature range would 
give little advantage as the damping would make accurate measurements impossible 
at still lower temperatures. ; 

There is only evidence of a positive dispersion in castor oil (Figure 4), but at 
rather lower temperatures than reported by Zachoval who found an increase of 1%, 
at 18°c. when increasing the frequency from about 2,000 to 6,000 kc/s., and the 


I 
1600 I 


» 1500 


1500 1300 


~25 0 25 ae 
Temperature °r ~85 -75 -65 =a 
4 Pee neers -50 -25 0 +20 Temperature °C 
Figure 4. Sound velocity in Temperature ‘c. 
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experimental 
——-— theoretical 


absolute value of the velocity at this temperature is somewhat higher than 
Zachoval’s, viz. 1,529 against 1,496 m/sec. at 2,000 kc/s. 

A negative disperison of the sound velocity was found in iso-butyl alcohol 
(Figure 5) over the whole temperature range. 


S45 .DISCUSSTON 


A negative temperature coefficient, at constant frequency, of the ultrasonic 
velocity of the type observed at higher temperature in all cases, and in iso-butyl 
alcohol at all temperatures within the experimental region, conforms with the 
results generally obtained with liquids other than water (Freyer, Hubbard and 
Andrews 1929, Prosorov and Nozdrev 1939, Suryaprakasan 1940, Willard 1947). 

The negative dispersion found at these temperatures, as shown in Figures 2, 3 
and 5, is small but it is nevertheless well substantiated, as pointed out above. This 
effect was quite unexpected and appears very difficult to account for. The only 
reference to the possibility of a negative dispersion is a theoretical paper by 
Ginsburg (1942) who bases his arguments on certain modifications of the Navier- 
Stokes equations, but the order of magnitude of the effect as predicted by this 
theory is much smaller than the values observed here. 


f 
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There is no singularity in the curves at temperatures and frequencies where 
dielectric relaxation occurs. The following table gives the temperatures in. 
degrees centigrade at which the loss factor has a maximum: 


Frequency (ke/s.) 2943 2102 689 
Geraniol * — 39° — 43° —51° 
Di-dihydrocitronellyl ether * —66° —72° 
Iso-butyl alcohol + —65° —81° 


* Schallamach 1946. 7+ Kobeko, Kuvshinsky and Shishkin 1938. 


It will be seen that there is no obvious connection between the two processes. 

The sharp rise of the sound velocity at lower temperatures appears to begin near 
the temperature where the curves cross and the positive dispersion starts. A. 
similar increase has been observed in liquid oxygen below —200°c. by Liepman. 
(1936) who, however, worked only at one frequency. In the case of our liquids 
it may well be that this effect is connected with the transition to a glass referred to. 
in §1. One would naturally expect a much higher sound velocity in the glassy 
state but it is surprising that the transition should make itself felt at comparatively 
high temperatures because a liquid like di-dihydrocitronellyl ether is brittle only 
in the neighbourhood of —140°c. Before discussing the possibility of frequency 
dependent mechanical constants, which is further suggested by the observed 
positive dispersion, the implications of the classical theory must be examined. 
The condition for a noticeable positive dispersion as predicted by hydrodynamics. 
is fulfilled here because the quotient (w7/)? can no longer be neglected in com- 
parison with unity. An analysis of the data is possible for the ether as the viscos- 
ities are known at the temperatures in question (Schallamach 1946). A certain 
complication arises from the negative dispersion which makes the difference 
between the velocities at the two frequencies ambiguous for this purpose. The 
relative differences are plotted against temperature in Figure 6, the shaded area 
around this curve indicates the experimental uncertainty, and the broken line gives. 
the theoretically expected differences calculated from equation (4), neglecting the 
very small correction at the lower frequency. It will be seen that there is rough 
agreement between the two curves but it appears that the experimental curve begins. 
to rise at a rather higher temperature than the theoretical curve, and that it falls 
below it at lower temperatures. This latter effect is more pronounced in castor oil 
for which the correction has been calculated with a viscosity extrapolated from data 
in the International Critical Tables. At —20-8°c., the observed relative velocity 
change is 1:5°% and the theoretical difference 4:6%. 

The experimental results acquire a different aspect it the liquid is taken to be a 
Maxwellian body characterized by a mechanical relaxation time 7 and a complex. 
viscosity 
gk =ng(l—jor)(Lte%). sea (6) 
Such a liquid has rigidity against alternating stresses the period of which is. 
comparable with the relaxation time. Introduction of a complex viscosity into 
equation (2) leads to expression (7) for the sound velocity which is correct to the: 
same degree of approximation as equation (4): 


2 (ow \?| « 1 1—w?r? ] 7 
pmol 3 (2) | eae tea aE Oe (7) 


where j1=7)/7 is the rigidity at infinitely high frequency. The second term in. 
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square brackets corresponds to the viscosity dispersion, which is now smaller than 
in an ideal liquid, but the first term gives the increase of the velocity due to the 
rigidity. It is difficult to estimate the value of the limiting rigidity » but the order 
of magnitude of the ratio «/u is approximately 3 for ordinary solids, and we feel 
justified in assuming that ps is, in any case, considerably smaller thanx. Equation 
(7) shows that the effective rigidity produces a greater dispersion than the viscosity. 
If, therefore, the conceptions underlying the formulation of the complex viscosity 
are accepted, it follows that any observed dispersion must be primarily due to the 
rigidity of the liquid. , 

The experimental evidence is insufficient to determine whether the dispersion 
reported here may be asscribed to complex viscosity. Measurements of the 
.absorption coefficient would most probably give an answer to this question. '' 

It may be argued that it is not necessary to restrict the discussion to a frequency- 
-dependent viscosity but that similar considerations could be applied to the com- 
pressibility. It must, however,.be borne in mind that this would involve the 
introduction of a compressional or ‘‘second”’ viscosity the existence of which is not 
recognized in hydrodynamics and which has by no means been proved yet 
experimentally. 
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REVIEWS OF BOOKS 


High Vacua, by Swami JNANANANDA. Pp. xiii+310. (New York: D. Van 
Nostrand Co. Inc., 1947.) 30s. 


It is a welcome feature of this book that prominence is given to the Kinetic Gas Theory, 
without a sure understanding of which it is impossible to understand many aspects of high 
vacuum technique. The author has, in fact, devoted one quarter of the book to this most 
important subject, which is treated under the following sectional headings : Introduction, 
Free Path Phenomena, Transport Phenomena, The Equation of State, Free Molecular 
Behaviour of Gases. The next quarter of the book is given up to a description of vacuum 
pumps, amongst them being “Pistonless High-Vacuum Pumps with Oils or Paraffins or 
High-Boiling Phlegmatic Liquids as Pumping Media ”’, i.e. Oil Diffusion Pumps. Then 
follows a long chapter giving descriptions of many types of gauges for pressure measurement, 
and the remaining three short chapters are devoted to the Technique of High Vacuum, 
Preparatory Operations for Vacuum Work, and finally, Production of High Vacua by 
Physical-chemical Methods. ‘The book is completed by the insertion of three short tables 
and an index. 

In appraising a book it is natural to ask why the book was written and what purpose it 
was intended to serve. ‘The answer to these questions is given in the last paragraph of the 
Foreword by George A. Lindsay: “The work supplies a real need for collected and 
organized information on the technique of the production of high vacua and should be a 
valuable reference for the worker in that field’’. It-is true to say, particularly with regard to 
the second or third chapters on the production and measurement of low pressures. 
respectively, which together form over half the book, that the author has carried out an 
immense amount of solid work in presenting an exhaustive catalogue of pumps and gauges, 
and from this point of view the aim of supplying a “‘ reference work ’’ has been achieved... 
But it is also pertinent to ask whether at this stage the worker in the high vacuum field requires 
a historical reference work or rather an up-to-date manual of vacuum practice, for it is quite 
certain that if a worker is looking for the latter he will not find it in this book. It is perhaps. 
the reviewer’s chief criticism that, whereas during the war there has been a phenomenal 
advance in the science of high vacuum technique, there is no reference in the book to any- 
work published later than 1940, and for this very reason the reviewer experienced an acute 
feeling of disappointment, which tended to overshadow the appreciation of the author’s 
work in compiling a history of the subject. This feeling of disappointment was unfortu- 
nately in no way lessened when a perusal was made of the manner in which the author dealt 
with certain aspects which are treated in most text-books on the subject. The reviewer has 
often found that in earlier works the explanation of certain features has been far from lucid, 
and, if one may be permitted to say so, this tradition has been well maintained in the present 
book. 

For example, the essential difference between Gaede’s “rate of exhaustion’, and 
Langmuir’s “‘ speed of a pump ”’ does not appear to be brought out at all in the explanation 
given on pp. 75, 76 where S has to do duty for both in equations (2-2 a) and (2-3 6). It is 
not very helpful to define the speed, as on p. 74, as the volume or mass of gas abstracted per 
unit time at a given pressure. 

Again, the discussion of conductivity and the flow of gases through pipes is extremely 
confusing. On pp. 77, 78 F is used to denote “‘ conductivity” or rate of flow, and Q the 
volume of gas flowing out, measured at a pressure of 1 dyne/em®. This should be compared 
with the definition on p. 103 for the output Q, viz. the ‘‘ product of the mean pressure and 
the volume removed per second, under a pressure difference of 1 dyne/cm?”’. In this case 
Q is the conductivity which was earlier denoted by F. On p. 228, Q is defined as “ the 
quantity of gas flowing per sec.”’ This is, of course, the best explanation, assuming that it is 
understood that ‘“‘ quantity ’’ refers to ‘‘ pressure X volume ”’, as indicated on p. 59. On 
p. 79, however, Q is defined as “ the quantity of gas taken out by the pump”. The 

‘confusion of thought is also apparent on p. 64 when F is incorrectly described as the rate of 
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flow in cm?/sec. per unit pressure difference, and again on p. 78 in the caption to figure 2-1. 
In transforming equation (1-129 a) to (1-129 6) the author has taken into account end effects 
only for the case of molecular flow, and has ignored them for the viscous flow. A com- 
parison of equations (1-126 b), which is correct, and (2-6), shows the presence of an alge- 
braical error, and this affects equations (1-129 b) and (1-129 ¢) equally. 

The author follows the usual procedure when he says on p. 115 that the speed of Lang- 
‘muir’s pump was fifty times as great as that of Gaede’s pump inferring that 80 cm*/sec. was 
the maximum speed to be obtained from a diffusion pump, and omitting any reference to 
‘Gaede’s later pumps described by him in 1923. 

Every reader may not agree with some of the statements made in the text. For example, 
it does not appear necessary in order to measure the speeds of pumps over an appreciable 
pressure range to use any other method than that described on p. 75. In some cases, 
‘however, the method described on p. 76 may be used as an alternative. Again, on p. 111, it 
js stated that in a diffusion pump, the gas molecules diffuse into the stream of mercury 
‘vapour, “‘ as the partial pressure of the gas is higher than that of mercury vapour ’’, and on 
p. 112, that ‘“‘ the optimum diffusion of gas molecules takes place only when the width of 
the aperture is comparable with their mean free path ”’. 

It is not the universally accepted way of sealing demountable joints to coat them both 
with a thin film of grease and then join them gently but firmly, thus filling the space between 
the surfaces with grease (p. 250) ; further, graphite is not a constituent of ““ Apiezon Q”’. 

There appear to be very few printer’s errors, although one or two have been observed, 
e.g. on p. 66 equation (1-129 a) should have the term D?/2:394L ; on p. 76 equation (2-4) 
should read loge [(p1—Po)/(Py—Po)| ; on p. 203, the equivalent of 0-4 mil. should read 
0-001 cm., and a similar correction is required on the top line of p. 204 ; on p. 217, equation 
(3-22) should read dr+dcf(p). R. WITTY. 


Technical Optics, Volume 1, by L. C. Martin. Pp. vii+393. (London: 
Sir Isaac Pitman & Sons, Ltd., 1948.) 40s. net. 


In the present dearth of books on optics the appearance of the first volume of a revised, 
enlarged (and renamed) edition of Professor Martin’s Introduction to Applied Optics is 
more than welcome. The scope of this work, as that of its forbear, is wide—almost too _ 
wide for anything like completeness to be attempted, although this disadvantage is offset 
‘by the very generous list of references at the end of each chapter. 

The first four chapters treat the paraxial theory of lens systems, the wave theory of 
light (including an introduction to the diffraction theory of aberrations), and the aberrations 
of symmetrical systems. The later chapters provide self-contained introductions to the 
optics of vision, physical optics and doubly refracting crystals, the elements of the theory 
of radiation and dispersion, and finally chapters dealing with the working and testing of 
‘optical surfaces and with the theory of spectacles. 

In a work embracing such a wide field it is gratifying that a careful reading only revealed 
‘one error of consequence. In dealing with the theory of Petzval curvature (p. 137) the 


usual expression [: G = ¥) | is obtained, but it is stated that this result is not valid in 
the case when the object is at the centre of curvature of the refracting surface. The proof 
given earlier (p. 23) in support of this contention refers to the curvature of the surface 
‘containing the tangential focus and not to the Petzval surface. In fact the usual form of 
the expression for the Petzval curvature remains valid no matter what the position of the 
object. A proof valid in this special case is given in Whittaker’s Theory of Optical Instru- 
ments (Cambridge Mathematical Tracts), although the proof given there breaks down when 
the object and the refracting surface coincide. It is a surprising fact that no universally 
valid proof of this well-known theorem ever appears to have been published. 

Apart from this, one can have nothing but praise for this book. One can scarcely say 
the same for the price. This, following the general upward trend, will no doubt make 
many students think twice about buying a book which they ought to buy without any 


hesitation. H. H. HOPKINS. 
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ABSTRAGTS FOR SECTION A 


Collision Problems and the Theory of Radiation Damping, by J. HAMILTON. 


ABSTRACT. The relation between radiation damping theory and the problem of the 
radiationless collision of an electron with a field of force is discussed. It is shown that 
the divergence difficulties, such as appear in the damping treatment of Compton scattering, 
do not arise, and the “ cut-off’ of damping technique is not necessary. ‘The same result 
applies to the relativistic collision of two electrons. 


Damping Theory and the Propagation of Radiation, by J. HAMILTON. 


ABSTRACT. It is shown that a careful examination of the concepts involved in radiation 
damping theory makes it possible to treat the propagation of energy from one atom to 
another through the electromagnetic field. An exact and continuous solution of the 
propagation problem is given. 


Continuous y-Emission in Neutron- Proton Collisions, by M. Kroox. 


ABSTRACT. ‘The cross-section for Bremstrahlung in neutron—proton collisions is 
calculated with a central exchange-force model and for energies up to 20 Mev. Its value 
is found to lie between 10-28 and 10-?° cm?. This probably is too small to be measured 
with experimental techniques at present available. 


Resonance Transfer of Electronic Energy in Organic Crystals, by E. J. Bowen, 
E. Mixiewicz and F. W. Snirtu. 


ABSTRACT. ‘The resonance transfer of energy in crystals of organic compounds when 
excited by visible or ultra-violet light has been examined by measuring the fluorescence 
of the systems naphthalene—anthracene, anthracene—naphthacene, naphthacene—pentacene, 
anthracene—acridine, anthracene-phenazine, and anthracene—di- and hexabromobenzenes. 
‘The conditions for “‘exciton’”’ trapping and the relations of resonance transfer to fluores- 
cence quenching are discussed. 


Further Experiments with an Adjustable Geiger-Miuller Counter, by A. G. FENTON 
and E. W. FULLER. 


ABSTRACT. An adjustable Geiger-Miiller counter and special electronic circuits have 
been used to study the effect of several counter variables on the occurrence of multiple 
pulses. The results show that the plateau characteristics combine the effects of an increase 
in efficiency and an increase in the probability of a multiple pulse as the operating potential 
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israised. It was found that the probability of occurrence of multiple pulses does not depend 
simply upon the charge per pulse, but that it depends also upon the anode diameter and the 
total pressure of the filling. j 

Curves have been drawn showing the charge per pulse as a function of the operating 
potential. A change of slope occurs in these curves when the charge per pulse is approxi- 
mately equal to the charge on the active portion of the anode. 

The counter was used as an ionization chamber at operating potentials below the Geiger 
region. It was found that as the potential is raised above the threshold of the proportional 
region the ionization current varies with potential as log Joc V while at higher potentials the 
relation is adequately expressed by log log Ja V. Near the Geiger region a sharp break 
occurs in the log log curve, which is taken to indicate that the discharge begins to spread at 
this stage. 


The Steady Non-Uniform State for a Liquid, by R. EISENSCHITZ. 


ABSTRACT. The Smoluchowski equation for the distribution of the relative coordinates 
of a representative pair of molecules in a liquid is formulated in such a way that the effect 
of non-uniform velocity or temperature is taken into account. ‘The equation is reduced 
to an ordinary linear differential equation of the second order. One of the two constants 
of integration is readily determined whereas’ the determination of the other constant may 
lead to inconsistencies. "The molecules are assumed to interact with the average potential 
that is appropriate to thermal equilibrium. If it is assumed that this potential becomes 
infinite of the first order at distance zero then the problem of viscous flow and the problem 
of thermal conduction are uniquely soluble. Approximate evaluation of the transport 
coefficients shows the right kind of temperature effect. The result is in accordance with 
the previous treatment by means of the cell model. 


Dislocation Theory of Yielding and Strain Ageing of Iron, by A. H. CoTTrReLi 
and B. A. Bitsy. 


ABSTRACT. A theory of yielding and strain ageing of iron, based on the segretation of 
carbon atoms to form atmospheres round dislocations, is developed. The form of an 

atmosphere is discussed and the force needed to release a dislocation from its atmosphere is 

roughly estimated and found to be reasonable. ‘The dependence on temperature of the 

yield point is explained on the assumption that thermal fluctuations enable small dislocation - 
loops to break away; these loops subsequently extend and cause yielding to develop 

catastrophically by helping other dislocations to break away. The predicted form of the 

relation between yield point and temperature agrees closely with experiment. 

Strain ageing is interpreted as the migration of carbon atoms to free dislocations. The 
rate of ageing depends upon the concentration of carbon in solution and the estimated 
initial rate agrees with experiment on the assumption that about 0-003°% by weight of 
carbon is present in solution. 


Investigations on the Binding Energy of Heavy Nuclei, by R. Husy. 


ABSTRACT. The binding energy of a heavy nucleus is one of the quantities which a 
satisfactory nuclear force theory ought to predict correctly. Approximate calculations of 
this binding energy are performed, taking as nuclear force the static interaction of Moller 
and Rosenfeld. The parameters appearing in the interaction have been fixed by data on 
light nuclei. 

The Fermi-gas model is taken as basis for the calculations, a first solution being obtained 
by a first-order perturbation (or variation) method. For refinement, a second-order pertur- 
bation calculation is made, and Svartholm’s variation-iteration method is attempted. 
The first-order results yield about 10% of the required energy, and the second-order about 
40% ; a reasonable prediction of the size of the nucleus is obtained. The reliability of the 
perturbation results is uncertain ; Svartholm’s method does not appear to be well suited to 
the investigation of heavy nuclei. = 

Two alternative modifications are made to the nuclear force to make it yield the correct 
binding energy according to the perturbation calculation. 
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